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ABSTRACT 

In a series of papers, we aim at stepping towards characterizing physical properties of the AGN dust torus by combining IR high-spatial 
resolution observations with 3D clumpy torus models. In this first paper, we present mid-IR imaging and 8 - 13yum low-resolution 
spectroscopy of 9 type I and 10 type 2 AGN. The observations were carried out using the VLT/VISIR mid-IR imager and spectrograph 
and can be considered the largest currently available mid-infrared spectro-photometric data set of AGN at spatial resolution SIOOpc. 
These data resolve scales at which the emission from the dust torus dominates the overall flux, and emission from the host galaxy 
(e.g. star-formation) is resolved out in most cases. The silicate absorption features are moderately deep and emission features, if seen 
at all, are shallow. The strongest silicate emission feature in our sample shows some notable shift of the central wavelength from 
the expected 9.7 /jm (based on ISM extinction curves) to ~10.5/jm. We compare the observed mid-IR luminosities of our objects to 
AGN luminosity tracers (X-ray, optical and [O iii] luminosities) and found that the mid-IR radiation is emitted quite isotropically. 
In two cases, IC 5063 and MCG-3-34-64, we find evidence for extended dust emission in the narrow-line region. We confirm the 
correlation between observed silicate feature strength and Hydrogen column density which was recently found in Spitzer data at lower 
spatial resolution. In a further step, our 3D clumpy torus model has been used to interpret the data. We show that the strength of the 
silicate feature and the mid-IR spectral index a can be used to get reasonable constraints on the radial dust distribution of the torus 
and the average number of clouds A'o along an equatorial line-of-sight in clumpy torus models. The mid-IR spectral index a is almost 
exclusively determined by the radial dust distribution power-law index a and the silicate feature depth mostly depends on A'o and the 
torus inclination. A comparison of model predictions to our type I and type 2 AGN reveals that average parameters of a = -I.O + 0.5 
and A'o = 5-8 are typically seen in the presented sample, which means that the radial dust distribution is rather shallow. As a 
proof-of-concept of this method, we compared the model parameters derived from a and the silicate feature strength to more detailed 
studies of full IR SEDs and interferometry and found that the constraints on a and A'o are consistent. Finally, we might have found 
evidence that the radial structure of the torus changes from low to high AGN luminosities towards steeper dust distributions, and we 
discuss implications for the IR size-luminosity relation. 
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1. Introduction 

The dust torus is one of the key ingredients of the unifica tion 
scheme of AGN (lAntonuccilll99l lUrrv & Padovanilll995h . It 
must be pointed out, however, that in recent years the torus pic- 
ture has evolved away from a "donut" towards a more general 
circumnuclear, geometrically- and optically-thick dust distribu- 
tion. The dust reprocesses the optical/UV photons of the accre- 
tion disk and re-emits the received energy in the infrared (IR). 
Thus, the torus can be directly studied using the IR emission of 
AGN, unless the AGN is radio-loud. Within the last years, in- 
vestigations based on IR SEDs characterized the torus emission 
and confirmed the basic picture of the unification scheme. Most 
notable are the silicate features at 10 and 18//m which are gen- 
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erally seen in absorption in type 2 AGN and in emission in type 
1 objects. 

However, more detailed studies - mainly using data from 
the IRS and MIPS instruments on-board the Spitzer satel- 
lite - revealed subtle differences to model predictions. In par- 
ticular, silicate emission features in type 1 AGN are much 
weaker than expected from early torus models which used 
smooth dust distributions. As shown in the literature, this 
might be explained by dumpiness of the dust in the torus 
which leads to an overall suppression of the emission fea- 

2005; 



ture ([Nen kova et al. 2002; Dullemond & van Bemmel 
et al. 2006; Schartmann et al...2008: .Nenkova et alJ I 
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Honig & Kishimot^l^l^ However, as will be sho wn in the 
second paper of this series ("Honig & Kishimotol20ld , pre-print 
available on astro-ph), weak emission features are not a generic 
property of clumpy torus models but might be used to constrain 
some model parameters of the dust distribution in clumpy tori. 
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Table 1. Characteristics of our type 1 AGN mid-IR spectroscopic sample 



Object 


Type" 


z 


Ol" 


log Luminosities (erg/s) 






Object scales 










(Mpc) 




Ly 


^-[Ollll 




'■sub * (pc) 


8 yum (pc)'' 


res. (''sub) 


NGC 3227 


S1.5 


0.0049 


20.4 


42.40 (22.80) 


42.20 


40.51 


42.63 


0.014 


24.2 


1680 


NGC 3783 


S1.5 


0.0108 


44.7 


43.19 (22.47) 


42.98 


41.30 


43.57 


0.035 


52.5 


1482 


NGC 4593 


Sl.O 


0.0101 


42.0 


42.89 (20.30) 


42.88 


40.69 


43.18 


0.032 


49.5 


1560 


ESQ 323-G77 


SI. 2 


0.0159 


66.3 


42.76 (23.25) 


43.91 


41.05 


43.75 


0.104 


77.1 


741 


MCG-6-30--15 


S1.5 


0.0087 


35.8 


42.86 (21.67) 


43.41 


40.36 


43.09 


0.058 


42.3 


727 


IC 4329A 


S1.2 


0.0170 


70.5 


43.87 (21.65) 


44.11 


41.21 


44.20 


0.131 


82.0 


625 


MARK 509 


S1.5 


0.0335 


141 


43.94 (20.70) 


44.01 


42.29 


44.25 


0.116 


158 


1365 


NGC 7213 


S1.5L 


0.0051 


21.2 


42.33 (20.60) 


42.76 


39.87 


42.50 


0.028 


25.2 


913 


NGC 7469 


S1.5 


0.0151 


62.7 


43.33 (20.61) 


43.09 


41.44 


43.87 


0.040 


73.1 


1808 



2007) 



200S ) 



— Notes: " AGN types from I Veron-Cettv & VeronI i2006h : S: Seyfert, L: LINER; * Luminosity distance based on CMB referenc e fram e 
redshifts from NED a nd Hq = 73, Q.^ = 0.27, and flvac = 0.73; Lx = absorption-corrected 2-lOkeV luminosity from iDadina 
except for NGC 3227 jTueller et al.ll200§) . logarithmic intrinsic Hydrogen column densities A'h in units of cm"^ from lTueller et al 
except ESO 323-G77 (average of 'Beckmann et al.' (2006) and' Malizia et all ( l2007h ); Ly = vL ,,(5 500 A); from'Be ntz et al.l l l2009h . except 
NGC 7213 (Winkler 1992), ESO 323-G77 (Schmid et al. 2003), and MCG-6-30-15 dBennert e t al. 2006). Additional extinction correction 
apphed to IC 4329A and ESQ 323-G77 I Wink ler et alj|1992n . and NGC 3783 (Ay = 1.0 f rom Ha/H6; . de Griio et al. 1992); " [Oin] total 
line fluxes from Schmitt et a ll r2003) (NGC 3783 NGC 4593, MCG-6-30-15, NGC 7213). lMelendezet al . (2008) (NGC 3227, IC 4329A, 
MARK 509, NGC 7469), and Schmid et al.|( l2003h (ESO 323-G77); ' Li2^„, = vL,,(12yum); this work; sfrom the fit to /C-band reverberation 
mapping data as shown in iKishimoto el all ( 120071) ; observed U - J or U - K time lag radii: 0.071 pc (NGC3783), 0.084 pc (MARK509), 
0.044 pc (NGC7469), Spatial scale corresponding to the theoretical diffraction limit for a 8.2 m VLT telescope: 0'.'25 at 8//m, 0'.'31 at 10//m, 
a'37 at 12yum; 



One of the problems of Spitzer is the comparably low spa- 
tial resolution of ~3" at lOyum. On the other hand, the torus 
near- and mid-IR emission originates from the inner few parsecs 
around the AGN which corresponds to < (X'l even for the near- 
est Seyfert galaxies. As a result, the Spitzer data is prone to con- 
tamination from the host galaxy or circumnuclear star-formation 
which can be identified by strong PAH emission lines. One way 
to overcome this problem is a decomposition of the data into sev- 
eral components. However, this requires good knowledge about 
generic SEDs of each component, which is, of course, difficult 
when aiming at characterizing the torus emission, and introduces 
additional parameters which leads to parameter degeneracies. 

The most direct way of studying the dust torus is IR inter- 
ferometry. Recently, the IR emission source in a small num- 
ber of nearby type 1 and type 2 AGN have been resolved in 
the m id-IR using the interferometric instrument MIDI at the 
VLTI ('Jaff^e et al. 2004; Tristram et al. 2007; Beckert et al."2008^ 
iRabaneta l. 2009; XristrameLak 2009; Burtscher et al. 2009i) 
and in the near-IR using the K eck interferometer (Swa in et al.l 
l2003t iKLshimoto. el al.l l2009ah . While interferometry provides 
the most direct ac cess to the toru s structure and brightness dis- 
tribution (Kishimoto e l al.ir2009a ). observations are limited to 
the brightest AGN which are in reach of current facilities. A 
compromise between Spitzer and interferometry can be achieved 
by using ground-based 8 m-class single telescopes. Such obser- 
vations usually do not have the potential of firmly resolving 
the torus, but may isolate the nuclear dust e mission from any 
contaminating source in the host g alaxy (e.g. iHorst et al.ll2009l: 
iGandhi et aI.l2009HMason et al.l20 09'). Following this approach, 
we observed 19 nearby AGN using the mid-infrared imager and 
spectrograph VISIR at the ESO Very Large Telescope (VLT) 
Paranal Observatory. In this second paper of our series, we 
present results from mid-IR spectro-photometric observations of 
our nearby AGN sample and interpret the observations with our 
3D c lumpy torus model (iHonig et al.ll2006l; iHon ig & Ki shimotol 
12010 ). Details about the type 1 and type 2 sub-samples will be 
presented in Sect. |2] The observations and data reduction are 
described in Sect. [3] Sect. |4] discusses the results of the obser- 



vations. In Sect.|5] we analyze the mid-IR characteristics of our 
sample and compare it to AGN properties. The data is then in- 
terpreted using our 3D clumpy torus model in Sect. |6] We sum- 
marize our main conclusions in Sect. |7] 

2. Properties of the AGN sample 

2.1. Sample selection 

Our main goal is to obtain the highest spatial resolution spectro- 
photometric data set of AGN yet obtained in the mid-IR. 
However, the sample cannot be considered "complete" in any 
respect, and restrictions and selection criteria are outlined in the 
following. However, we consider the sample to be "typical" (or 
representative) for the Seyfert galaxy population in the Galactic 
vicinity (see below). Moreover, it should give us an idea of the 
"clean" nuclear emission at scales of tens of parsecs. 

The original idea for selecting objects was based on a de- 
mand of high-spatial resolution. In particular, we aimed for ob- 
jects which have a resolution <100pc around 10 /vm using the 
8.2 m UT3-telescope at Paranal. Thus, we are limited to AGN 
at angular-diameter distances <70 MpcQ. Since the ESO/Paranal 
observatory hosting VISIR is located in the Southern hemi- 
sphere, the whole sample is limited to objects mostly at Southern 
declinations. In addition, for easy execution of the observa- 
tions in service mode, the AGN were selected to be brighter 
than 100 mJy in most of the A^-band. Based on these crit eria, 
we mined the AGN catalog of IVeron-Cettv & VeronI (l2006l) and 
compared the objects to previous low-spatial resolution mid-IR 
data from Spitzer and ISO to be able to estimate which ob- 
jects are bright enough for VISIR. Since our ultimate goal con- 
cerns typical nearby AGN, i.e. Seyfert galaxies, all peculiar (e.g. 
LINERs) were skipped. To avoid synchrotron contamination, we 
also excluded radio-loud objects. Finally, since we are interested 
in the characteristics of the AGN, we avoided any objects where 

' Please note that the distances given in Table [T] & |2] are luminosity 
distances D^. The corresponding to our 70 Mpc limit would be < 
73 Mpc. 
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Table 2. Characteristics of our type 2 AGN mid-IR spectroscopic sample 
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Object 


Type" 


z 




log Luminosities (erg 


ys) 




Object scales 










(Mpc) 


Lx (A'h) 


-'^[OlllJ 




'■sub * (pc) 


8 yum (pc)^ 


res. (fsub) 


NGC2110 


Sli 


0.0080 


33.2 


42.60 (22.57) 


40.35 


43.00 


0.028 


39.2 


1392 


ESO 428-G14 


S2 


0.0063 


26.0 


. . . (>24) 


40.62 


42.73 


0.020 


30.8 


1509 


MCG-5-23-16 


Sli 


0.0095 


39.5 


43.24 (22.47) 


40.63 


43.53 


0.052 


46.6 


900 


NGC 4507 


Slh 


0.0128 


53.0 


43.34 (23.46) 


41.57 


43.69 


0.062 


62.1 


1004 


MCG-3-34-64 


Slh 


0.0176 


73.2 


42.78 (23.69) 


41.45 


44.14 


0.104 


84.9 


814 


NGC 5643 


S2 


0.0047 


19.2 


. . . (>24) 


40.55 


42.49 


0.016 


22.9 


1471 


NGC 5995 


S1.9 


0.0257 


107 


43.52 (21.95) 


42.01 


44.06 


0.095 


122 


1294 


IC 5063 


Slh 


0.0109 


45.3 


42.87 (23.28) 


41.28 


43.82 


0.072 


53.2 


737 


NGC 7582 


Sli 


0.0044 


18.3 


42.04 (22.98) 


40.15 


42.71 


0.020 


21.8 


1093 


NGC 7674 


Slh 


0.0277 


116 


43.72 (24.00) 


41.93 


44.26 


0.120 


132 


1103 



— Notes: " AGN types from lVeron-Cettv & VeronI i2006h : S: Seyfert, L: LINER, Sli: optical type 2 with broad Hydrogen lines in the infrared; 
Slh: type 2 with detected polarized optical broad lines; * Luminosity distance based on C MB reference f rame redshifts fro m NED and 
Ho = 73, = 0.27, and flyac = 0.73; Lx = ab sorption-corrected 2-lOkeV luminosity from Dadinal d2007'). except IC 5063 ( Tuel ler et al.l 
|2008), MCG-3-34-64 and NGC 4507 jHonig e t al. 2008a), NG C 5995 dH orst et al. 2008), and NGC 7674 (Bianchi et al. 2005); logarithmic 
intrinsic Hydrogen column densities in units of 10^- cm"^ from Tueller et al. (2008), except ESO 428-G14 ('Delia Ceca et al. 2008), NGC 5643 
{Maiolino et al. 1998), and NGC 7674 (Bianchi et al. 2005, changing-look AGN); [Ora] fluxes from Schmittetal. (2003) (NGC 4507, 
IC 5063, NGC 7674), Melendez et al. (2008) (NGC 5643, NGC 7582). iGu et'al.i ( 12006 ^ (ESO 428-G14, MCG-3-34-64), Haas et al. (2003) 
(NGC 21 10), .Wink ler(1992) (MCG-5-23-16), and Tran (200 3) (NGC 5995); " Liz^m = vL,(l2 ^um); this work; f based on eq.[D; « Theoretical 
diffraction limit for a 8.2 m VLT telescope: 0'.'25 at 8 jim, Qf.'3 1 at 10 yum, 0'.'37 at 12 yum; 



the nucleus is heavily obscured by host-galactic dust lanes, e.g. 
as in Circinus. Thus, any obscuration pattern seen in the data of 
most of our objects (i.e. silicate absorption features) is suppos- 
edly intrinsic to the nuclear environment. However, we acknowl- 
edge that IC 4329A and NGC 7582 have dust lanes passing over 
the nucleus, so that part of the observed absorption properties 
may originate in the host galaxy (in particular for the Seyfert 2 
NGC 7582; see Sect. |6]for possible consequences on our analy- 
sis and modeling). In summary, we selected 9 type 1 and 10 type 
2 AGN. 

2.2. Observed and intrinsic scales 

In Table [T] & |2] we list the basic characteristics of our type 1 
and type 2 AGN sub-sample, respectively. In case the objects are 
Compton-thin, the intrinsic (absorption-corrected) 2- lOkeV X- 
ray luminosity can serve as a proxy for the total luminosity of the 
AGN. For the type 1 AGN (see Table[TJ, we also provide optical 
luminosities which are more direct tracers for the accretion-disk 
luminosity. Observed scales for each objects are provided at the 
reference wavelength of 8 fim. 

Another way to look at our main requirement of high-spatial 
resolution is not the observed scale but the intrinsic scale of 
each object. Since most of the nuclear mid-IR radiation in radio- 
quiet Seyfert is presumably coming from dust emission, the fun- 
damental scaling relation between the dust sublimation radius 
and the AGN luminosity, rjub L}^^, provides such a distance- 
independ ent scaling. If w e take UV/IR-reverberation mapping 
data (e.g. lSueanuma et an i2006. and references ther ein) and use 
the calibrated relation from IKishimoto el al.l (l2007h . we obtain 

Tsub = 0.36 pc X ( Lv / 10 erg/s j . For the type 1 AGN, we can 
directly calculate the expected rsub from their optical luminos- 
ity and compare rjub to the spatial resolution of our observations 
(see Table[T]i. As can be seen, all of the objects which meet our 
resolution selection < lOOpc (corresponding to Di < 73 Mpc) 
have an intrinsic spatial resolution < 2 000 rsub- Since the intrin- 
sic scale is the ultimate factor determining how good the AGN 
can be isolated from the host galaxy, it is probably safe to in- 
clude additional objects in the sample which have a spatial reso- 



lution >100pc, but still an intrinsic resolution < 2 000 rsub with- 
out compromising our goal of highest spatial resolution. Thus, 
we include Markarian 509 which adds to the higher luminosity 
end of the type 1 sub-sample. 

Since the UV/optical emission from the accretion disk is 
thought to be obscured in type 2 AGN, a direct comparison be- 
tween observed scales and intrinsic scales is not possible via 
the reverberation-based size-luminosity relation. However, us- 
ing our type 1 sub-sample, we are able to convert the size- 
luminosity relation from the optical to other wavebands. Since 
we are dealing with mid-IR observations, the most convenient 
way is a conversion of the size-luminosity relation to 12//m. 
This wavelength is mostly outside the silicate feature so that it is 
not affected by possible anisotropics due to silicate absorption or 
emission. Recent high-spatial resolution studies (also using the 
VISIR instrument on the 8.2-m VLT/UT3 telescope) have shown 
that type 1 and type 2 AGN follow basically the same Lmir -^-a'- 
relation in the luminosity range which is covered by our sample 
(iHorst et al.ll2008l; iGandhi et al.ll2009l) . Anisotropy between the 
two samples are within the observational errors and are smaller 
than about a factor of 2 to 3. Thus, Li2/jm can also serve as a 
proxy for the AGN luminosity of our objects and enables us to 
estimate intrinsic scales. We note that this is only a good approx- 
imation as long as high spatial resolution data is available. 

As a first step, we determine the correlation between Ly 
and L[2f,m in the type 1 sub-sample, assuming that the cov- 
ering factor is similar for all objects. We find that logL^ = 
(6.9 + 9.6) + (0.84 + 0.22) x logLn^,™ (Spearman rank 0.88, null- 
hypothesis probability 5 x 10"^). Within errors, this is consistent 
with Lv L\2f,m, with a ratio Lv/Lnfim = 0.59^q34, which we 
assume in the following. Using this correlation, we obtain a scal- 
ing relation for our AGN sample of 

.,.. = (0.2rs:.>cx(^)"' (1, 

Despite the uncertainty of a factor of 1.5, this should enable us 
to give at least an estimate of intrinsic scales of our type 2 ob- 
jects. Intrinsic scales based on eq. ^ are given in Table |2] As 
can be seen, while NGC 5995 and NGC 7674 have observed 
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Table 3. VISIR mid-IR photometry of our 9 type 1 AGN measured by Gaussian fitting to science target and calibrators. In case 
several calibrators were available, the fluxes and FWHM are mean values. 



Object 


Filter 


/I, 


AA 


Flux 


FWHM 


Observing date 




name 


(yum) 


(//m) 


(miy) 


target/calib 




NGC 3227 


Arlll 


8.99 


0.14 


179.6 ± 


11.3 


a'32/a'31 


2008-03-22X02:37 




PAH2ref2 


11.88 


0.37 


320.1 ± 


21.9 


a'37/a'33 


2008-03-22X02:41 


NGC 3783 


PAHl 


8.59 


0.42 


323.9 ± 


33.7 


a'3i/a'26 


2008-03-20X06:08 




Arlll 


8.99 


0.14 


350.9 ± 


23.2 


a'30/0'.'29 


2008-03-20X06:16 




SIV 


10.49 


0.16 


558.4 ± 


23.2 


a'3i/a'3i 


2005-04-17X01:27 




PAH2ref2 


11.88 


0.37 


595.8 ± 


39.2 


0'.'34/a'32 


2008-03-20X06:20 




Nellrefl 


12.27 


0.18 


685.0 ± 


37.2 


a'35/0'.'32 


2005-04-17X01:36 


NGC 4593 


Arlll 


8.99 


0.14 


181.8 ± 


13.8 


0'.'31/0'.'26 


2008-04-01X06:32 




PAH2ref2 


11.88 


0.37 


278.8 ± 


20.4 


a'36/a'3i 


2008-04-01X06:35 


ESQ 323-G77 


Arlll 


8.99 


0.14 


298.9 ± 


11.5 


a'40/a'29 


2009-05-19X00:58 




PAH2ref2 


11.88 


0.37 


377.4 + 7.6 


a'37/a'36 


2009-05-10X01:01 


MCG-6-30-15 


SIV 


10.49 


0.16 


311.7 ± 


11.4 


a'32/0'.'38 


2006-04-14X03:58 




PAH2 


11.25 


0.59 


359.6 ± 


24.7 


a'33/0'.'33 


2006-04-14X04:03 




Nellrefl 


12.27 


0.18 


377.5 ± 


11.2 


0'.'36/0'.'36 


2006-04-14X04:07 


IC 4329A 


Arlll 


8.99 


0.14 


743.4 ± 


22.0 


a'30/a'29 


2009-05-19X01:12 




PAH2ref2 


11.88 


0.37 


1014.1 d 


: 18.8 


a'36/0'.'36 


2009-05-10X01:15 


MARK 509 


SIV 


10.49 


0.19 


207.3 ± 


17.0 


a'35/0'.'34 


2006-06-14X10:20 




PAH2 


11.25 


0.59 


240.8 ± 


27.1 


a'32/a'32 


2006-06-14X10:25 




Nell 


12.81 


0.21 


243.0 ± 


19.8 


0'.'30/0'.'35 


2006-06-14X10:30 


NGC 7213 


SIV 


10.49 


0.16 


239.1 ± 


22.0 


a'33/0'.'32 


2006-07-14X09:58 




PAH2 


11.25 


0.59 


273.3 ± 


34.5 


a'37/0'.'34 


2006-07-14X10:03 




Nellrefl 


12.27 


0.18 


245.6 + 


16.9 


a'35/a'32 


2006-07-14X10:07 


NGC 7469 


SIV 


10.49 


0.16 


447.8 ± 


15.6 


a'39/a.'3i 


2006-07-12X07:21 




PAH2 


11.25 


0.59 


470.1 ± 


38.8 


a'42/a'34 


2006-07-12X07:25 




PAH2ref2 


11.88 


0.37 


505.9 ± 


25.4 


a'34/a'32 


2006-07-14X09:29 




Nellrefl 


12.27 


0.18 


595.0 ± 


18.0 


a'36/a'34 


2006-07-12X07:30 




NeIIref2 


13.04 


0.22 


629.7 ± 


16.6 


a'35/a'47 


2008-08-26X06:48 



scales >100pc, their intrinsic scales are < 2 000 r^ab, even when 
accounting for the factor of 1.5 uncertainty in eq. ([TJ. These two 
objects have been added to our sample to have higher luminosity 
Seyferts in both the type 1 and type 2 sub-samples. We note that 
the relation based on Lx instead of Li2fim would not improve the 

results (r,„b = 0.45!0:32pc x (Lx/10« evg/sf'\ 

In summary, our objects can be described as a sample of typi- 
cal nearby Seyfert galaxies - obscured and non-obscured - span- 
ning the luminosity range from about 10"^^ erg /s to 10'^* erg/ s. 
The intrinsic spatial resolution achieved by our observations is 
< 2 000 rsub- Please note that we accidentally included one op- 
tically broad-line (BL) LINER in our sample, NGC 7213 (al- 
legedly a "type 1 LINER"), which will be treated as a type 1 
AGN in the rest of the paper but discussed separately where ap- 
plicable. 



3. Observations and data reduction 

3.1. Spectroscopy 

We used the VISIR mid-infrared imager and spectrograph 
mounted on the 8.2 m UT3 telescope at the ESO/Paranal obser- 
vatory in Chile. The observations have been carried out in ser- 
vice mode in ESQ periods 78, 80, 82, and 83. In total, 9 Seyfert 
I and 10 Seyfert 2 galaxies have been observed in low-spectral 
(LR) resolution mode (R ~ 300). The standard configuration for 
LR long-slit spectroscopy results in a pixel resolution of 0'.'127, 
which samples the PSF in the A^-band by ~2.5 pixels. To cover 
the full A^-band, 4 different spectral settings have to be used with 
central wavelengths at 8.5, 9.8, 11.4, and 12.4//m. We used a 



slit width of 0'.'75 which is approximately 2-3 times larger than 
the FWHM achieved with VISIR (a.'25 - a.'39 across the A^- 
band) and minimizes the risk of slit losses. All of our targets are 
unresolved point-sources, except for NGC 7469 where some off- 
nuclear emission is expected from the circumnuclear star-burst 
ring. 

For the initial steps of data reduction, we followed the 
standard ESO pipeline. The individual chopped nodding cycles 
(chopping width 8", nodding in parallel mode) were combined 
and wavelength-calibrated using the Common Pipeline Library 
(CPL) recipes in ESOREX. After that, we used our own pro- 
cedures to calibrate and extract the spectra. First, the remain- 
ing sky-offset has been removed. For that, we fitted a low-order 
polynomial in spatial direction to each wavelength bin (typically 
of orders to 2). We note that this procedure can only be used 
because we are not interested in very large-scaled smooth struc- 
ture — which would not be visible anyway because of the flux 
limit. This polynomial sky subtraction flattened the image back- 
ground significantly (see Fig.[T]i. 

However, a periodic background pattern in spatial direction 
remained in each wavelength bin, with a frequency length of 16 
pixels and not depending on chopping width, frequency, or posi- 
tion angle. In usual observing conditions, this pattern dominates 
the total variance of the background. The ratio of object peak to 
background variation peak can be as low as 4: 1 in AGN, thus be- 
ing very disturbing when dealing with faint mid-IR objects. We 
developed an efficient method to remove this background pattern 
using a "periodic background map (PBM)". The PBM is gener- 
ated by creating a cube with several copies of the sky-removed 
science array, each copy shifted by 16 pixels in spatial direction 
with respect to the previous copy. Finally, all shifted science- 
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Table 4. VISIR mid-IR photometry of our 10 type 2 AGN measured by Gaussian fitting to science target and calibrators. In case 
several calibrators were available, the fluxes and FWHM are mean values. 



Object 


Filter 


A, 


AA 


Flux 


FWHM 


Observing date 




name 


(yum) 




(miy) 


target/calib 




NGC2110 


PAHl 


8.59 


0.42 


156.6 + 9.9 


a' 33/0'.' 32 


2006-10-10X07:28 




Arlll 


8.99 


0.14 


169.1 ± 6.4 


a'31/a'29 


2006-10-10X07:31 




PAH2ref2 


11.88 


0.37 


294.3 + 8.0 


a'33/0'.'35 


2006-10-10X07:40 


ES0 428-G14 


PAHl 


8.59 


0.42 


112.7 ± 18.5 


a'37/0'.'41 


2007-03-01X03:27 




Arlll 


8.99 


0.14 


122.8 ± 14.4 


a'4o/a'3o 


2007-03-01X03:34 




PAH2ref2 


11.88 


0.37 


207.1 ±9.4 


a'43/a'32 


2007-03-01X03:42 


MCG-3-34-64 


SIV 


10.49 


0.16 


591.9 ±22.4 


a'39/a'25 


2006-04-09X05:23 




PAH2 


11.25 


0.59 


685.8 ±51.0 


a'43/0'.'31 


2006-04-09X05:31 




Nellrefl 


12.27 


0.18 


889.3 ±21.2 


a'44/a'34 


2006-04-09X05:35 




NeIIref2 


13.04 


0.22 


1116.0 ± 60.7 


a'46/0'.'36 


2006-01-18X07:58 


NGC 4507 


SIV 


10.49 


0.16 


489.2 ±43.3 


a'32/0'.'32 


2006-04-15X02:58 




PAH2 


11.25 


0.59 


554.7 ±43.6 


a'32/a'33 


2006-04-15X03:02 




Nellrefl 


12.27 


0.18 


631.4 ± 10.8 


a'34/0'.'34 


2006-04-15X03:07 


MCG-5-23-16 


PAHl 


8.59 


0.42 


337.2 ± 10.5 


a'27/0'.'33 


2007-01-30X06:09 




Arm 


8.99 


0.14 


358.2 ± 16.7 


a'28/a'32 


2007-01-30X06:32 




PAH2ref2 


11.88 


0.37 


633.4 ± 24.0 


a'33/a'35 


2007-01-30X06:40 


NGC 5643 


Arlll 


8.99 


0.14 


143.3 ± 33.5 


a'35/a'27 


2009-05-03X04:58 




PAH2ref2 


11.88 


0.37 


287.1 ±9.1 


a'37/0'.'35 


2009-05-03X05:01 


NGC 5995 


SIV 


10.49 


0.16 


266.0 ± 10.7 


a'37/0'.'38 


2006-04-14X04:48 




PAH2 


11.25 


0.59 


332.5 ± 17.7 


a'38/0'.'31 


2006-04-14X04:52 




Nell 


12.81 


0.21 


381.8 ± 20.2 


0'.' 40/0'.' 33 


2006-04-14X04:57 


IC 5063 


SIV 


10.49 


0.16 


608.6 ± 22.0 


a'47/a'37 


2006-05-05X09:54 




PAH2 


11.25 


0.59 


726.6 ± 24.7 


a'43/0'.'38 


2006-05-05X09:59 




PAH2ref2 


11.88 


0.37 


924.9 ±25.0 


a'42/0'.'38 


2005-06-10X08:06 




Nellrefl 


12.27 


0.18 


1036.5 ± 56.2 


a'43/0'.'37 


2006-05-05X10:03 


NGC 7582 


PAHl 


8.59 


0.42 


327.2 ± 37.7 


a'51/0'.'41 


2008-08-18X09:29 




Arm 


8.99 


0.14 


235.9 ±27.2 


a'33/a'29 


2008-08-23X08:41 




PAH2ref2 


11.88 


0.37 


383.6 ±26.2 


a'38/a'33 


2008-08-23X08:38 


NGC7674 


Nell 


12.81 


0.21 


518.4±21.7 


a'45/a'38 


2006-07-13X08:55 



array copies are combined by median -filtering each pixel. The 
resulting PBM contains only the 16-pixel-frequency background 
without flux from the (point-like or slightly resolved) science 
target. In Fig. [T] we illustrate the PBM removal for the 12.4 yum 
setting of NGC 7213. The solid light-gray line represents the raw 
data after sky subtraction (mean over all rows). For these data, 
the PBM has been determined, as shown by the dotted dark-gray 
line. The dashed black line shows the final data after PBM re- 
moval. As can be seen, the noise variation is significantly sup- 
pressed as compared to the original data. 

The described procedure has been applied to both science 
and calibrator data. After removing all background, wavelength- 
dependent conversion factors have been determined from the 
calibrators and the science data were flux-calibrated accordingly. 
We refrained from airmass corrections since the differential air- 
masses between target and calibrator were rather small, so that 
corrections would be within the calibration errors. In the case 
of MCG-3-34-64, we found a beam-centering problem. In fact, 
visual inspection of acquisition images revealed that only part of 
the object was placed within the Qf.'lS slit. We used the acquisi- 
tion images in the N_SW filter (made through the open slit) of 
the science target and the calibrator to quantify the loss due to 
this problem. The measured flux for MCG-3-34-64 in this filter 
is F(N_SW) = 42 1.06+ 24.75 mJy (wavelength: 8.85 + 1.35 yum). 
Comparing to the integrated flux in the spectrum, we found that 
only about 70 + 4% of the object flux was located within the slit. 
Thus, we corrected the flux-calibrated spectra of MCG-3-34-64 
accordingly. 
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Fig. 1. Illustration of the periodic background map (PBM) re- 
moval for the 12.4;um setting of NGC 7213. The solid light- 
gray line represents the raw data after sky subtraction (mean over 
all rows). The PBM is shown by the dotted dark-gray line. The 
dashed black line presents the final data after PBM removal. 
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Fig. 2. VISIR low-resolution spectroscopy (black solid line and gray error bars) and photometry (red filled circles) of 8 type 1 AGN. 
Overplotted (blue dashed line) are Spitzer data with approximately 10 times less spatial resolution. The hatched areas mark regions 
with strong sky lines which are difficult to calibrate. Prominent mid-IR emission line positions are indicated. The slit position angle 
(RA^) is given in the lower right corner of each panel. 




3.2. Imaging photometry 

Along with the spectroscopy data, we acquired VISIR images 
in several mid-infrared filters to compare and validate the abso- 
lute and relative flux calibration of our spectra. Together with 
archival VISIR data, all of our objects were observed in at least 
2 different narrow-band filters between 8 and 13/im, except 
NGC 7674 (only 1 filter). For photometry, the standard CPL 
pipeline products were used for both science target and cali- 
brator The imaging was performed in parallel or perpendicular 
chop/nod mode, so that 3 or 4 beams appear on the detector for 
the science data (parallel: 2 negative single-beams and 1 posi- 
tive double-beam). On the other hand, the standard stars were 
observed almost exclusively in perpendicular chop/nod mode re- 
sulting in 4 single-beams (2 negative, 2 positive beams). 

Each science and calibrator beam was fitted by a 2- 
dimensional Gaussian to resemble the core of the Airy disk of 
the PSF. We confirmed that our science targets are point sources 
in the mid-IR by comparing the Gaussian width of the tar- 
gets with the standard stars. The measured Gaussian FWHM of 
science target and associated calibratoj^ are given in Table [3j 
Although the science targets usually have slightly larger FWHM 
than the calibrators, the differences are not significant, except for 
MCG-3-34-64 and IC 5063 which we will discuss in Sect.O 
in detail. In most cases, the differences between the FWHM of 
the science targets and calibrators are probably the result of PSF 
instabilities of VI SIR which are known to be a common problem 
dHorst et al.l 12009). Therefore we consider that our sources are 
mostly unresolved. After the Gaussians have been removed, all 
science targets showed the first Airy ring without any additional 
emission source, except for NGC 7469. This galaxy has a well- 
known nuclear star-burst ring within 2" of the nucleus which we 
spatially resolve in our observations and which contains more 
flux than the first Airy ring. 

The final fluxes were obtained by calculating conversion fac- 
tors from the integrated intensity of the Gaussian fits of the 
calibrators and using these factors for the respective integrated 
Gaussian intensity of the science targets. Since we are only in- 

^ Where several calibrators were available, the given calibrator 
FWHM is the average of all calibrators. 
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Fig. 3. VISIR low-resolution spectroscopy (black solid line and 
gray error bars) and photometry (red filled circles) of the LINER 
galaxy NGC 7213. Overplotted (blue dashed line) are Spitzer 
data with approximately 10 times less spatial resolution. The 
hatched areas mark regions with strong sky lines which are dif- 
ficult to calibrate. Prominent mid-IR emission line positions are 
indicated. The slit position angle (PA) is given in the lower right 
corner. 



terested in point source fluxes, this method of optimized extrac- 
tion has the advantage that any flux from off-nuclear emission is 
excluded. This provides the best possible estimate for the torus 
emission which is presumably unresolved. The resulting fluxes 
are presented in Tables |3] & S All data was obtained with a 
chop throw of 8", except for the archival PAH2ref2 filter data 
of NGC 7469 and the archival NeIIref2 filter data of MCG-3- 
34-64 where a chop throw of 10" was used. 
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Fig. 4. VISIR low-resolution spectroscopy (black solid line and gray eiTorbars) and photometry (red filled circles) of 10 type 2 AGN. 
Overplotted (blue dashed line) are Spitzer data with approximately 10 times less spatial resolution. The hatched areas mark regions 
with strong sky lines which are difficult to calibrate. Prominent mid-IR emission line positions are indicated. The slit position angle 
(PA) is given in the lower right corner of each panel. 
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Fig. 4. — continued. 



4. Results and discussion 

4. 1 . VISIR spectra and photometry 

In Figs. |2]&[3]we show the individual VISIR 8-13 /urn spec- 
tra of our type 1 AGN (8 Seyfert galaxies and 1 BL LINER) 
together with results from VISIR photometry. In the lower-right 
corner of each panel, we noted the slit position angle (PA) for 
each object. The corresponding data of the 10 Seyfert 2 galaxies 
are shown in Fig. |4] In addition to these data, the correspond- 
ing Spitzer IRS low-resolution spectra are shown for the same 
wavelength range. The spatial resolution of the IRS spectra is 
~ 3" while the VISIR data resolves scales of ~ (X^ at lO/zm. 
For NGC 7582, no Spitzer spectrum was available yet. The po- 
sitions of prominent emission lines seen in AGN are shown. The 
gray-hatched areas in each spectrum mark regions of extensive 
sky-line emission/absorption which are difficult to calibrate due 
to some degree of variability within the night and/or compro- 
mised S/N ratios. 



All VISIR spectra (except MCG-5-23-16; see below) have 
lower or equal fluxes than the Spitzer spectra, which is expected 
for the higher spatial resolution. A similar result for VISIR 
photometry on a sampl e of type 1 and type 2 AGN has re- 
cently been reported by iHorst et aP |2009). Together with the 
fact that the imaging-photometry is well consistent with the 
spectro-photometry, this confirms self-consistency of the ex- 
tracted fluxes and supports consistency of our calibration ap- 
proach with respect to other observations. 

In MCG-5-23-16 the VISIR spectrum is slightly above both 
the Spitzer spectrum and the VISIR photometry data point. Since 
the VISIR photometry matches the Spitzer fluxes we account 
this discrepancy to some unknown calibration error in spec- 
troscopy (e.g. a slightly deviating spectral shape of the cali- 
brator). In the extreme case around 12jum, the difference is 
50mJy or 7% of the flux. In NGC 7469, the VISIR photome- 
try is slightly but systematically lower the VISIR spectrum. The 
reason for this is that the extraction window used for the spec- 
troscopic data includes some small degree of emission from the 
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Fig. 5. Spitzer IRS minus VISIR differential spectra of NGC 3227 (left), NGC 5643 (middle), and NGC 7469 (right). The differential 
spectra are compared to the scaled-down IRS spectrum of M82 (blue-dashed line) which is commonly used as a template for 
extragalactic star-formation. 



surrounding star-burst ring (see also Sect. 14.21 and iHorst et al] 
(2009) ) while our optimal flux extraction for imaging isolated 
the nuclear point source. 

As most ground-based mid-lR instruments, VlSlR uses 
chopping and nodding to subtract the strong sky emission. 
However, this procedure does not only affect the background but 
also has an influence on the science data: If an emission region 
is extended on scales larger than the chop-throw, this emission is 
strongly reduced. The degree of suppression depends on the ac- 
tual extension and spatial flux gradient. It is strongest for objects 
which are more or less uniformly emitting light. On the other 
hand, point-like objects or emission regions confined within the 
chop-throw are not affected. Thus, chopping can be considered 
a useful tool to reduce or suppress extended host galaxy emis- 
sion. Such contribution from the host galaxy can be stellar light 
or emission from extended star-formation projected onto the nu- 
cleus, in particular if the host galaxy is strongly inclined. In con- 
clusion, our spectra and photometry can be expected free of ex- 
tended host galaxy emission, containing only emission from the 
point-like AGN and its immediate vicinity. 

4.2. PAH emission and star-formation activity around the 
AGN 

Thirteen of the Spitzer IRS spectra of our objects show more 
or less pronounced PAH emission features at 11.3 fim, and in 
some cases the 7.7/8.6 //m PAH complex is also well visible. 
Since PAH emission is associated with star-formation, the pres- 
ence of these features in the IRS data indicate that star-formation 
is occurring in the central ~ 3" of the AGN or, due to the lack 
of chopping, is projected onto the nucleus in the Spitzer data. 
Common to all of our high-spatial resolution VISIR spectra as 
compared to the Spitzer IRS data is the lack, or at least strong 
suppression, of PAH line emission. Most objects do not display 
any 11.3/im PAH emission, while in NGC 3227, NGC 5643, 
NGC 5995, NGC 7469, and NGC 7582 some minor PAH fea- 
ture can be seen. For Mrk 509 the situation is inconclusive due to 
the low quality of the 12.4yum setting, and for NGC 4593 we do 
not have data in the corresponding spectral setting. In the cases 
with some little 11.3yum PAH emission, it is difficult to judge 
if some marginally remaining feature at around 8.6 yum is also 
present (e.g. by comparing to the Spitzer data), given that this 
feature is generally broader in wavelength and weaker It has to 
be pointed out that the suppression of the PAH emission features 



is not an effect of different spectral resolution but only caused 
by the different spatial resolution and observing techniques. 

In Fig.|5] we show examples of IRS minus VISIR differential 
spectra of objects with prominent PAH emission in the IRS spec- 
tra, NGC 3227, NGC 5643 and NGC 7469. For that, the VISIR 
data were downgraded to the same spectral resolution as IRS 
{R ~ 100). The differential spectra are compared to the scaled- 
down version of the IRS spectrum of M82 (blue-dashed line), 
which is often used as an extragalactic star-formation template. 
The match for all three sources is reasonably good, especially for 
NGC 7469. Thus, we conclude that the differential spectra from 
spatial regions between ~ 0'.'3 - 3" are predominantly showing 
star-formation emission in the continuum and lines. 

The main differences of spectra taken with Spitzer and 
VISIR are spatial resolution and observing technique. Since our 
spatial resolution is about a factor of 10 better than the IRS data, 
the emission regions producing the PAH features can simply be 
located at scales between 0'.'3 and 3". On the other hand, only 
in few cases individual star-forming regions are actually seen 
in the VISIR ima ges, and image s of these cases have recently 
been shown by Ho rst et al.l (l2009l) . The most notable example is 
NGC 7469 with its well-known star-burst ring at about 2" from 
the nucleus. When integrating over a 0'.'75 x 3" aperture, we re- 
cover part of the IRS spectrum. The remaining "missing flux" 
probably originates from extended emission outside the VISIR 
aperture or from the host galaxy projected onto the nucleus. 
Since VISIR uses the chopping/nodding technique, any extended 
host galaxy emission at scales beyond the chop throw will be 
eliminated from the data or, at least, significantly reduced. As 
a result, the nuclear point source flux is free of contaminating 
emission from the host galaxy, even if part of the host emission 
falls onto the nucleus. 

Still, most of the PAH emission in Spitzer is supposed to 
come from the vicinity of the AGN. Since we selected galaxies 
with only moderate inclination, projection effects should play a 
minor role. Consequently, we see a significant reduction of the 
PAH emission from scales of 1 kpc down to <100pc. What is the 
reason for the suppression of the PAH emission features? PAH 
dust grains are prone to photo-destruction by high-energetic pho- 
tons as emitted by an AGN or in its vicinity (e.g. from hot ther- 
mal plasma). Therefore it is quite plausible that they become 
less abundant at smaller distances from the AGN. On the other 
hand, a lack of PAH emission can also point to reduced or no 
star-formation activity. While it is not possible to distinguish be- 
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tween both effects only from the PAH emission, the reduction 
of PAH emission features from IRS to VISIR goes along with 
a reduction in continuum flux in our spectra. If the only rea- 
son for suppressed PAH emission were photo-destruction of the 
associated grains, we would expect that the continuum level re- 
mains roughly constant. Thus, we conclude that the dominating 
effect of suppressed PAH emission is an actual decrease in star- 
formation activity at smaller distances from the AGN. 



4.3. The N-band silicate feature in AGN at h/gh spatiai 
resoiution 

Silicate absorption and emission features are the most evident 
spectral signatures in the A^-band. One of the surprising discov- 
eries of Spitzer was based on observations of silicate emission 
features in type 1 AGN. The features have been much weaker 
than expected from early toru s modeling (e. g. Pier & Kroli 



2.0 



1993 



iGranato & Danesg Il994t lEfstathiou & Rowan-Robinso 
19951) . The same also holds for silicate features in absorp- 
tion as seen in type 2 AGN, although very deep features 
are sometimes seen in ULIRGs or type 2 AGN where the 
host galaxy presum ably contributes by a large fraction to 
the obscuration (e.g. iLevenson et aT]l2007t IPoiletta et al.ll2008t 
iMartinez-Sansigre et al.l 20091) . Since the Spitzer data has a spa- 
tial resolution of several arcseconds in the A^-band, it was not 
clear initially to what degree the weakness of the silicate fea- 
tures was a resolution effect. Recent ground-based mid-IR spec- 
troscopy of single objects at high spatial resolution, however, 
suggested that even at resolution better than 1" silicate ab- 
sorption and e mission features were moderate (e.g. iRoche et"al] 
l2007tlMason et al. 2009) . Here we present a much larger sample 
of type 1 and type 2 AGN observed at sub-arcsecond resolution 
so that the silicate feature characteristics can be studied more 
systematically (see Sect. 15.3b . 

Fig.|2]illustrates that most type 1 AGN exhibit a rather weak 
emission feature, if we can detect any at all. A small bump in 
the continuum can be seen somewhere between 9 and 12yum, 
e.g. in NGC 3783 and MCG-6-30-15. The other features are 
only revealed when plotting log vFy or fitting the continuum (see 
Sect. 15.3b and some seem to be rather absorption than emis- 
sion features (e.g. NGC 3227). The strongest silicate emission 
features is displayed by the weak or LINER AGN NGC 7213 
(Fig. |3]l. Based on the spatial resolution of our data of <100pc 
in all objects except MARK 509, we conclude that the weakness 
of the silicate emission features in type 1 AGN is not a spatial 
resolution effect but an intrinsic property of the features and pos- 
sibly of the emission region. 

The type 2 AGN show a variety of feature characteristics (see 
Fig.lDi. While most type 2s have moderate silicate features in ab- 
sorption, the NGC 4507 A^-band spectrum resembles a feature- 
less type 1 spectrum. In NGC 2110, we even see the silicate fea- 
ture in emission. It is not surprising seeing such different char- 
acteristics in type 2s since we expect various torus inclinations 
reaUzed in these objects, ranging from moderate to very high ob- 
scuration ("edge-gracing" to edge-on geometries). This is illus- 
trated by the variety of Hydrogen column densities A^h observed 
in X-rays for these AGN (see Table |2|. In fact, NGC 2110 has 
the lowest A^h within the type 2 su b-sample and infrared broad 
emi ssion lines have be en reported dVeron-Cettv & Veronll2006l) 
(see lMason et al.ll2009l for possible scenarios in NGC 21 10). 

One of the mysteries of the silicate feature in AGN is an 
actual or apparent shift of the emission feature towards longer 
wavelengths. Absorption features in type 2 AGN are centered 
very close to 9.7 fiva, just as expected from opacity curves of 
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Fig. 6. Continuum-normalized silicate emission features of 
PG1211-H143 (red) and NGC 7213 (blue). A linear continuum 
has been fitted to both objects based on their 8.3 and 12.7 jum 
fluxes and the observed spectrum was divided by this continuum 
fit. The small inset shows both features in log plot normalized 
for their different feature strength. As a comparison we overplot 
the extinction coefficient of standard ISM dust with Ossenkopf 
silicates (black-dotted line). The blue-dashed and red-dotted- 
dashed lines denote the central wavelengths of NGC 7213 and 
PG 1 2 1 1 + 1 43 respectively. 



the Galactic ISM dust (IChiar & TielerisI l2006h . Clumpy torus 
models based on ISM dust well reproduce the overall IR dust 
re-e mission and the sili c ate absorption features in detail (e.g. 
Hon ig et all 120061 120071: iPoUetta et all l2008t ISchartmann etall 
20081) . For type 1 AGN, the situation is less clear. Early re- 
ports of the silicate emission feature in a number of quasars 
discussed that the central wavelength might be shifted towards 
longer wavelength s (Hao et al. 20 05; Sigb enmorgen et al.|l2005J; 
ISfiirm et alj2005h . Recendy, Nik utta et all (2009) suggested that 
radiative transfer effects might cause the silicate emission fea- 
ture in AGN to shift towards longer wavelengths. In their anal- 
ysis, they require a certain number of clouds to intervene the 
line-of-sight of other clouds to cause some absorption tip at the 
center of the silicat e feature. To occur at th e right wavelength, 
they suggest using lOssenkopf et alj d 19921) silicates with peak 
wavelength at lO.O/im instead of standard ISM. However, they 
note that any slightly larger number of clouds or different opac- 
ity would cause a noticeable absorption dip to be present in the 
center of the si l icate e mis sion feature, as shown in the models of 
iNenkovaetaP (12008') and'Honi g & Kishimotol (l2010 l). Based on 
our own modeling, we would expect that if such radiative trans- 
fer effects are the only reason for a shift or the silicate emission 
feature, then we would expect that analyzing a sample of objects 
would reveal both shifted features and those showing noticeable 
absorption dips within the emission feature. However, all of the 
type 1 AGN observed here and in addition NGC 21 10 as a type 
2 AGN with a silicate emission feature do not show such an ab- 
sorption dip within the emission feature. This is consistent with a 
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detailed study of the sili cate emission features o f 23 PG quasars 
observed with Spitzer (Schweitzer et al.' '2008*). The problem 
was also mentioned by Mason et al. (2009) when attempting to 
model the silicate emission feature in NGC 2110. In addition, 
we note that our modeling of silicate features with Ossenkopf et 
al. silicate opacities show both emission and absorption f eature s 
centered at lO.Oyum, in agreement with Nenkova e t all (l2008h . 
Since, however, the observed absorption feature c entral wave- 
length is at 9.1 um (see also iMason et a l. 2006; Roche et al.l 
l2007t iHonig et alJ l2007h . type 2 AGN would favor different 
opacity curves. 

In Fig. |6] we show a comparison between the continuum- 
normalized A^-band Spitzer spectrum of PG121 1 + 143 analyzed 
by Nikutta et al. (2009) and our VISIR spectrum of NGC 7213. 
To estimate the continuum, we made a linear fit to the 8.3 and 
12.7jum fluxes in Fy. The observed spectra are then divided by 
the continuum fit. Using this kind of fit ( which is to some ex- 
tent similar to the method suggested by ISirockv et al.l (l2008h 
and used by .Nikutta et al., (2009 )). we see that the peak emis- 
sion of PG1211-H143 is located indeed at lO.Ojum, while it is 
shifted to ~10.5 pvn. in NGC 7213. We compare these features to 
the extinction curve of standard ISM dust with Ossenkopf sili- 
cates and find that the PG121 1 + 143 silicate emission profile is 
well matched while a shift is evident in NGC 7213. Note that if 
we used the Spitzer IRS spectrum of NGC 7213 instead of the 
VISIR spectrum, the result would be more or less the same. 

In summary, we would conclude that simple radiative trans- 
fer effects due to absorption within the torus alone, either clumpy 
or smooth, are not capable of explaining the details of the sili- 
cate features and their characteristics in type 1 and type 2 AGN. 
There may be other effects contributing to the actual shape of 
the feature. One possibility is that there is some (radial) change 
of the grain size and/or dust composition based on the fact that 
large grains, preferably graphite grains, have a much higher sub- 
limation temperature than smaller grains, in particular silicates. 
Then, hotter regions might be dominated by graphite dust and 
larger grains while cooler regions have an ISM-like size and 
graphite-silicate mix, which may go along with a temperature- 
dependent c entral wavelen gth of the silicate feature (e.g. as il- 
lustrated in iKrugell l2008l Fig. 9.7). Consequently, exact de- 
tails on absorption and emission within the silicate feature are 
more complex than covered by recent torus models, although the 
essence of the features (i.e. feature strength) are well captured by 
these models. 



5. AGN continuum and dust emission properties in 
the mid-IR at <100 pc 

As shown in the previous section, the high-spatial resolution of 
VISIR allows us to isolate the nuclear mid-IR emission in AGN 
without too much disturbance from host-galactic sources. Thus, 
analysis of this data will reveal characteristics of the circumnu- 
clear region around the AGN. In the following, we will discuss 
several properties of the mid-IR continuum emission and the sil- 
icate features. NaiTow forbidden atomic emission lines which 
are also present in the spectra, like rAriii], [Siv ], and [Neii] 
have been already studied in part in iHonig et"an (^008a), and 
the whole sample will be dealt with separately in an upcoming 
paper. 
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Fig. 7. Absorption-corrected X-ray luminosities Lx of our sam- 
ple of AGN plotted against the 12 pm luminosity Li2;jm derived 
from the VISIR spectra, excluding Compton-thick objects (see 
Tables[T]&|2]i. The blue triangles are type 1 objects while the red 
squares are type 2 AGN. The dotted line is the best-fit correla- 
tion for our whole sample, while the dashed line represents a fit 
to the sample excluding the outliers (see Sect. 15. lb The gray sym- 
bols show the full sample of the latest version o f the Lmir - Lx- 
coiTelation as presented in lGandhi et al.l (l2009h . 



5. 7. Mid-IR continuum emission and its relation to AGN 
luminosity tracers 

The mid-IR emission is believed to originate from dust re- 
processing of the central accretion disk radiation. In particu- 
lar, optical/UV photons hit the dust in the torus and heat up 
the grains. The absorbed energy is then re-emitted in the IR 
and forms the general IR bump seen in multi-wavelength SEDs 
(e.g. lElvis et"al]ll99 4'). This redistribution of energy should pro- 
duce a fairly significant correlation between the X-ray (accre- 
tion disk tracer) and mid-IR (dust emission tracer) luminosity 
in AGN. Such a Lmir - Lx-correlation has been found by sev- 
eral auth ors using both space-based and grou i id-based observa- 
tions (e.g. Krabbe et al.'200lVLutz et al.l2004t lHorst et al.l2008l; 
lOandhi et al. 2009). One remarkable property of this correlation 
is the apparent isotropy among the AGN: Within the observa- 
tional scatter of the latest version of the coiTelation of about a 
factor of 2 to 3 (Gandhi et al. 2009), there seems to be no diff'er- 
ence between type 1 and type 2 AGN. 

In Fig.Qwe show the absorption-corrected X-ray luminosi- 
ties Lx of our sample of AGN plotted against the 12 pm lumi- 
nosity Li2^ni derived from the VISIR spectra (see Tables [T]&|2]|. 
Compton-thick objects are not shown in the plot. The errors of 
the mid-IR luminosities are <0.1 dex (typical photometric ac- 
curacy is much better than 15%). The X-ray luminosities are 
taken from single-epoch data, so that we consider typical un- 
ceitainties in Lx due to intrinsic variability of ~0.5 dex (see er- 
ror of the Gandhi et al. (2009) sample overplotted as gray sym- 
bols in Fig. |7]). Within the scatter of our observations, there is 
no apparent difference between type 1 and type 2 AGN. This 
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Fig. 8. [O III] luminosities of our sample compared to the 12yum mid-IR (left) and 2 - 10 keV X-ray luminosities. The blue triangles 
are type 1 AGN, the red squares denote type 2 objects. The dotted lines are best-fit correlations for the whole sample (left: L[Oni] 

r0.96±0.I2. „;„ut. r „ r0.85±O.I8\ 

, right. L[0 III] ocL^ ). 



is con siste nt with the larger sa mples analyzed by iHorst et aP 
(120081) and lGandhi et al.l (120091) : the latter sample is shown for 
comparison as gray symbols in Fig. |7] Note that some ob- 
jects in the original Gandhi et al. sample are duplicated in our 
sample. However, different X-ray luminosities have been used 
and the mid-lR fluxes were extracted from our VISIR spec- 
tra. Using our much smaller sample leads to a nominal corre- 
lation Li2f,m Lx' °^''" '^ (Spearman rank pspearman = 0.82, null- 
hypothesis probabUity_4.92<10"^) which is consistent with the 
findings in lGandhi et al.l (1200 9*). There seems to be a small off'set 
of our sample towards lower mid-IR luminosities with respect to 
the Gandhi et al. study. This can be explained to a small degree 
by the fact that Gandhi et al. used 12.3 fim as the reference wave- 
length where fluxes are usually slighly higher than at 12.0/im 
as used here. Most of the difference, however, comes from the 
method used to extract the fluxes: Our p hotometry wa s mea- 
sured using the optimal extraction while in iGandhi eran(l2009l) 
aperture photometry was used which potentially includes some 
off'-nuclear emission. 

The non-detection of a difference between type 1 and type 
2 AGN within the error of observations might indicate a rel- 
atively large degree of isotropy of the AGN radiation in both 
the mid-IR and the X-ray. However, it has to be noted that the 
Lmir - Lx-correlation omits all Compton-thick objects. Unless 
one has a good idea of the intrinsic X-ray luminosity of these 
highly obscured obje cts (e.g. from Une emission at optically 
thin wavelengths as in! Honig et al.l2008ah . the observed isotropy 
must be considered as a lower limit. One such possibility is the 
commonly used [Om](/i5007 A) luminosi ty, although some de- 
gree of anisotropy has been reported (e.g. Netzer et ani2006h . In 
Fig. [8] we compare the [Oiii] luminosities of our sample to the 
12yum mid-IR and 2 - lOkeV X-ray luminosities, respectively. 
Both properties correlate well with L[Oni], but the ii2/jm-i'[0in]- 
relation is stronger by eye and by a statistical analysis. For our 
sample we find log L[o III] = (0.11+5.3)-i-(0.94+0.12)xlogLi2^m 



(Spearman rank p = 0.87, null-hypothesis probability 1.3x10 ^) 
and log L[o m] = (4.38 + 7.8) + (0.85 + 0. 1 8) x log Lx (Spearman 
rank p = 0.78, null-hypothesis probability 1.9 x 10""*). Note 
that for the correlation with Li2/jm, there are also Compton-thick 
objects included (NGC 5643 and ESO 428-G14). Thus, if the 
[Oiii] radiation is considered as being a good isotropy tracer, 
then Li2fim may be considered as being emitted quite isotropi- 
cally as well - at least more isotropically than X-ray radiation. 
This might also be interpreted as a sign that the mid-IR opti- 
cal depth of any obscuring medium (e.g. as traced by the X-ray 
hydrogen column density A^h) must be more transparent in the 
mid-IR, or that the mid-IR emission is emitted within or outside 
the X-ray-opaque medium. 

These findings seem to have a natural interpretation in the 
unification scheme: The dusty torus is both the X-ray-obscuring 
and the mid-IR-emitting medium. Recent studies have shown 
that the observed mid-IR isotropy (i.e. small dispersion between 
type 1 and type 2 AGN in the Lmir - L^-correlation) can be ex- 
plained within the fram ework of a clurnpy torus (e.g.fa onig et jTI 
2006; Horst et al. 2008: Nen kova et al. '2008VGa ndhi et al.l2009F 
Levenson et al.,,20 09). There is, however, an interesting point to 
consider: While IC 5063, ESO 323-G77, and MCG-3-34-64 
are outlier s in the Lmir - Lx -correlation (see Fig. Q and also 
iHorst et al . 2008; Gandhi et al. 2009), they are very close to the 
ii2//m-i'[0iii]-fit (see Fig.|8]l. This could be caused by absorption 
since all three objects have moderate A^h (logA^H ~ 23). On the 
other hand, all recent Lmir - L^-studies use A^n-corrected Lx so 
that it is difficult to imagine that the intrinsic Lx are still under- 
estimated by a factor of ~3-5. Another possibility would be that 
these objects are not "underluminous" in the X-rays but "overlu- 
minous" in the mid-IR, i.e. there is additional mid-IR emission 
other than the dust torus re-emission. This overhead emission 
is probably triggered by the AGN, since our spectra do not in- 
clude any major star-burst component (see Sect. 14.2b . Instead 
the additional mid-IR emission may originate from dust in the 
narrow-line region (NLR). The prime example for such a case 
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is NGC 1068. iMason et all (l2006h point out that the flux in the 
central I'.'l of this AGN (~80pc, which is a typical resolution 
element of our sample) contains only 30% contribution from the 
nuclear point source (presumably the torus) while 70% originate 
in extended emission, mostly from the NLR0. This coiTesponds 
to a mid-lR "overluminosity" of 0.3-0.4 dex in the Lmir - Lx- 
correlation, which is about the offset of the three sources. If we 
exclude these sources from the Lmir - Lx-correlation analysis, 
we obtain a tight relation of Li2^m <^ /^j^O-94±o.io (gpgarman rank 
Pspearman = 0.96, null-hypothesis probability 2.5 x 10"^), which 
is consistent within errors with Gandhi et al. (2009). 

An interesting aspect about this suggestion is that mid-lR 
VISIR imaging of the outliers MCG-3-34--64 an d IC 5063 actu- 
ally shows extended emission. |Horstetal] (l2009l) mention some 
"slight elongation" in their SIV, PAH2, and NelLl filter images 
of MCG-3-34-64, but attribute it to instrumental effects. We re- 
analyzed their images of MCG-3-34-64 and IC 5063, together 
with additional archival VlSlR images (see Fig.|9]l, by applying a 
2D Gaussian fit to both science target and calibrator The result- 
ing FWHM and position angles are presented in Table|5] MCG- 
3-34-64 shows at least 38-54% and IC 5063 shows at least 14- 
45% extension compared to the maximum elongation of the re- 
spective calibrators. All calibrators appear more or less round 
with axis ratios close to unity. On the other hand, the average axis 
ratios are 1.35 for MCG-3-34-64 and 1.26 for IC 5063. This 
goes in line with very similar position angles in each filter while 
the orientation of the axis in the standard star images is more ran- 
dom (average PA 54° for MCG-3-34-64 and 1 1 1° for IC 5063). 
In IC 5063, this is consistent with the orientation o f the extended 
[Oml emission as observed with HST (PA 115°: ISchmitt et af] 
l2003h . Despite no high-spatial resolution [Oiii] image is avail- 
able for MC G-3-34-64, we can compare our results to radio 
observations. ISchmitt et alj (12001 ') report linear extension of the 
nuclear radio emission (inner ~100pc) at 8.46 GHz towards PA 
39° which is only about 13° off from the mean elongation of 
the mid-IR emission. If we assume that the linear radio emission 
traces some jet-like emission and that the jet-like emission is ap- 
proximately in the same direction as the NLR, we can argue that 
the mid-IR emission in MCG-3-34-64 is also extended in the 
NLR direction. 

Actually, these 2 objects are the only ones in our sample 
where we see such a significant extension and elongation (see 
data in Tables [3] & 111 point-like images for some of our objects 
are shown in Horst et al. (2009)). Moreover, some members of 
our group attempt VLTI/MIDI mid-IR interferometry on both 
objects in the past. Although total fluxes were easily recorded, 
no fringes have been found. This is consistent with a very low 
12jum visibility (0.2-0.3 at about 60m baselines) in these ob- 
jects, meaning that 70-80% of the nuclear flux comes from ex- 
tended emission. Thus, based on (1) the absence of PAH emis- 
sion lines in our VISIR spectra, (2) the resolved emission regions 
observed in the mid-IR images of MCG-3-34-64 and IC 5063, 
(3) the significant elongation, (4) consistent position angles of 
the mid-IR extended emission region and the [O iii] emisison in 
IC 5063 and the linear radio emission in MCG-3-34-64, and (5) 
the low visibilities in mid-IR interferometry, we suggest a sce- 
nario where theses two objects have very mid-IR-bright NLRs - 
similar to NGC 1068 - which causes them to deviate from the 
Lmir - Lx-correlation. Whether or not this may also be an ex- 
planation for the apparent "overluminosity" o f AGN at the low- 
luminosity end of the Lmir - Lx -correlation dHorst et al.ll2008t 



' We note that the NGC 1068 data used in Figs /reffig:specindex to 
I12lare point source values from a C/.'4 extraction window. 
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Fig. 9. VISIR mid-IR contour images of MCG-3-34-64 (left 
column) and IC 5063 (right column). The different panels in 
each column represent images in different filters (SIV = forbid- 
den [Siv](10.51;um) line, PAH2 = 11.3//m PAH line, PAH2.2 
- 1 1.88 //m continuum, NelLl = 1 2. 27 //m continuum, NeII_2 = 
13.04yL/m continuum). Contour levels are scaled logarithmically 
from the peak in steps of 0.2dex down to the 4% peak flux level. 
Both objects show clear elongations, corresponding geometric 
properties based on Gaussian fits as listed in Table|5j 

iGandhi et al.l2009h has to be addressed by detailed mid-IR stud- 
ies. 

5.2. The mid-IR spectral index 

One fundamental observational property of the mid-IR emission 
is its spectral slope. Since a sizable part of the A^-band is affected 
by the silicate feature, the continuum slope has to be recovered 
from wavelengths outside the feature, which is centered around 
lOfim. Here we assume that fluxes at 8.5 /vm and 12.5 yum are 
mostly unaffected by the silicate feature (which we tested by 
inspecting ISM dust extinction curves and Spitzer IRS data cov- 
ering a larger wavelength range). 
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Fig. 10. Mid-IR spectral slope and spectral index for each object in our sample, plotted against the observed mid-IR luminosity 
Ln^im and the X-ray Hydrogen column density A^h- Type 1 AGN are shown as blue triangles, t ype 2 objects are rn arked by red 
squares. For reference, NGC 1068 is shown as an orange star, based on the mid-IR spectrum from lMason et al.l (l2006l) . 



Table 5. Gaussian fit geometric properties of MCG-3-34-64 
and IC 5063 base on VISIR images 



Filter 


object properties 


calibrator properties 




size 


PA 


size 


PA 




MCG-3-34-64 




SIV 


a'488 X a' 342 


59° 


a'316xa'298 


95° 


PAH2 




57° 


a'355 x0'.'320 


93° 


Nellrefl 


a'505 X a' 400 


45° 


^'365 xa'328 


82° 


NeIIref2 


a'562 X a' 435 


46° 


a'391 xa'384 


129° 




IC 5063 






SIV 


a'551 xa'417 


110° 


a'380xa'373 


175° 


PAH2 


a'508 X a'382 


109° 


a'379 xa'368 


69° 


PAH2ref2 


a' 452 X a' 394 


117° 


a'397 xa'385 


36° 


Nellrefl 


a'485 X a' 394 


107° 


a'384xO'.'380 


160° 



In Fig. [To] we show the spectral slopes (= flux ratios 
Fv(12.5//m)/Fv(8.5;um)) and spectral indices Fy oc y" (based 
on the flux ratios) for all of our objects. In the left panel, we 
plot the spectral slope log Fv(12.5/vm)/Fv(8.5yum) against the 
observed mid-IR luminosity. The distribution is rather wide, so 
that a clear correlation cannot be seen. However, there might be 
a slight trend that the reddest objects in both the type 1 and type 
2 sub-sample are on the lower luminosity end of the plot. What 
can clearly be seen is that type 1 and type 2 AGN do not show a 
huge difiference in the spectral index a, i.e. type Is can be as red 
as type 2 AGN, although a marginal difference may be detected. 
The nominal mean spectral indices are cfAGNi = -1.65 ± 0.44 
for type 1 AGN and OAcm = -2.07 + 0.54 for type 2 AGN. 
This almost similarity in spectral index becomes even more evi- 
dent when comparing the Hydrogen column density A^h with the 
spectral index. There is at best a marginal trend of the spectral 
index with obscuration towards the AGN - the most obscured 
objects appear to be redder than the mean in this sample. On the 
other hand, NGC 1068 as a Compton-thick object has a - -1.5 
in the mid-IR, so that it is placed in the bulk of the other objects 



Table 6. Properties of the VISIR mid-IR spectra. For each ob- 
ject we list the 8.5-12.5/im spectral index a, the silicate feature 
strength Fsi/Fc, and the corresponding optical depth in the sili- 
cate feature TsiUcate (see Sects. l572l & l5.3l for details) 



Object 


a 


logFsi/Fc 


^silicale 


NGC 1068 


-1.52 


-0.16 


0.38 


NGC 21 10 


-1.66 


0.01 


-0.03 


ESO 428-G14 


-2.62 


-0.25 


0.57 


MCG~5-23-16 


-2.36 


-0.14 


0.32 


NGC 3227 


-2.22 


-0.06 


0.13 


NGC 3783 


-1.84 


-0.03 


0.06 


NGC 4507 


-1.68 


0.02 


-0.05 


NGC 4593 


-1.51 


0.06 


-0.15 


ESO 323-G77 


-0.81 


0.02 


-0.05 


MCG-3-34-64 


-2.06 


-0.09 


0.20 


MCG-6-30-15 


-1.53 


0.03 


-0.08 


IC 4329A 


-1.46 


-0.01 


0.02 


NGC 5643 


-3.05 


-0.20 


0.45 


NGC 5995 


-1.34 


-0.05 


0.13 


MARK 509 


-1.24 


0.11 


-0.25 


IC 5063 


-2.38 


-0.12 


0.28 


NGC 7213 


-2.21 


0.10 


-0.23 


NGC 7582 


-1.46 


-0.45 


1.03 


NGC 7469 


-2.01 


0.01 


-0.02 


NGC 7674 


-2.11 


-0.08 


0.19 



(see star in Fig.[TO] left). The spectral index for each object has 
been tabulated in Table |6] 

The fact that type 1 and type 2 AGN have very similar spec- 
tral indices might seem quite surprising. The first-order assump- 
tion based on the torus picture would be that type 1 AGN are, on 
average, bluer than type 2 AGN, because most of the IR emis- 
sion we see in type 1 s comes from hot dus t in the inner part of 
the torus. On the other hand, as shown in iHonig & Kishimotol 
(2010, e.g. Sects. 3.2, 3.4 & 3.5.1, and Figs. 6, 7, 9, 10, 11 & 
12), the way the dust is distributed around the AGN can have a 
much stronger eff'ect on the mid-IR properties than inclination 
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Fig. 11. Silicate feature analysis of our AGN sample observed with VISIR. Type 1 AGN are shown as blue triangles, type 2 AGN 
are marked by red squares. Some objects have been i dentified by name. For reference, we include NGC 1068 with mid-lR properties 
based on the Gemini spectrum bv lMason et al.l 1*2006). Left: Dependence of the silicate feature strength Fsi/F^ (see text for details) 
on the mid-lR luminosity Li2;jm- The dashed line shows the nominal relation for type 2 AGN, the dotted line is the corresponding 
relation for type Is. Middle: F^i/Fc plotted against the spectral slope and spectral index of the mid-lR continuum emission. The 
dashed line shows the nominal relation for type 2 AGN, the dotted line is the corresponding relation for type Is. Right: F^-JFc 
and observed optical depth in the silicate feature TsiUcate shown as a function of the Hydrogen column density A^h- The dashed line 
represents the correlation F^-JFc oc 0.12 x logA^n- Gray circles are type 1 data points if the central wavelength of the silicate 
emission feature is 10.3 yum instead of 9.8 yum. The corresponding correlation Fsi/F^ oc 0. 14 x log A^h is shown as dotted line. 



effects - in particular if the torus is clumpy so that transitions 
in obscuration properties from type 1 to type 2 are smooth. We 
demonstrated that if the radial dust distribution is steep (rj,- oc /' 
with a - -1.5 ... - 2.0), i.e. most of the dust is confined to 
small radii, the resulting SEDs are quite blue in both typ e 1 an d 
type 2 orientations of the torus (see Honig & Kishimoto 20101 , 
Figs. 6 & 9). On the other hand, more shallow dust distributions 
(a = -0.5 ... - 1 .0) show redder colors for all orientations. Thus, 
in the framework of a clumpy torus, the similarity in mid-lR 
spectral indices can be naturally explained. A larger difference 
in spectral indices, however, may emerge at wavelengths shorter 
than 8 micron. 



5.3. The silicate feature and its relation to AGN properties 



In Sect 14.31 we discussed qualitatively the general appearance 
of silicate features at high spatial resolution. Here, we aim for 
a quantitative analysis of the silicate feature and its relation to 
other observed properties. For that, we fit a linear continuum to 
the 8.5 fim and 12.5 fim flux of each object. From that we inter- 
polate the expected linear continuum flux F^ at 9.8 /im around 
the center of the silicate feature (but see Sec. 14.31 for details on 
the central wavelength). Then, we take the ratio between ob- 
served flux and fitted continuum flux, Fsi/F^, as a measure for 
the stren gth of the silicate feature, which is similar to the strategy 
used by Hao et al] (l2007h . The silicate feature strength of each 
object is listed in Table |6] 

Fig. [TT] shows the silicate strength plotted against several 
observed parameters. The right axis shows the translation into 
an observed optical depth, Tsiucate, within the silicate feature. It 
has to be pointed out, however, that this does not reflect the 
actual optical depth along the line-of-sight towards the central 
engine, but is a result of the radiative transfer within the torus 
(see also lHonig & Kishimoto! I201Q) . In the left panel, we show 



the silicate strength as a function of mid-lR luminosity Li2fim- 
As in all plots, the type 1 AGN are marked by blue triangles, 
and red squares are used for type 2 AGN. For reference, we put 
NGC 1068 in all of th e plots based on the Gemini spectrum of 
the nucleus reported in lMason et al.l ( l2006l) . The type 1 AGN are 
all clustered around log Fs JFc = 0, reflecting the absence or ex- 
tremely shallow silicate features. They do not show any depen- 
dence of the silicate feature on Ln^im in the covered luminosity, 
with a nominal relation log Fsi/F^ <x (0.00 ± 0.03) x logLn^/m 
(dotted line in the left panel of Fig.fTTT). The type 2 AGN show 
a much larger variety of silicate feature strengths from moderate 
absorption to slight emission, as seen also in Fig.lD The deepest 
feature is observed in NGC 7582. This galaxy shows a well- 
known star- burst ring close to the n ucleus which is at high incli- 
nation (e.g. ' Wold & Gallianoll2006l) . It seems quite likely that a 
large fraction of the obscuration in this object is not intrinsic to 
the torus but caused by dust located in or around the star-burst 
ring or the inner host galaxy. In general the silicate absorption 
features seem to become shallower with luminosity in the type 
2 AGN. The dashed line in the left panel of Fig. [TT]is the nom- 
inal fit to the type 2 sub-sample, logF^JF^ oc (0.12 + 0.06) x 
logLi2^m- We would not call this a correlation (Spearman rank 
p - 0.58, null-hypothesis probability 8.2 x 10"^), but the data 
suggest at least that the silicate feature shows stronger differ- 
ences between type 1 and type 2 AGN at lower than at higher 
luminosities. If this is really true, it indicates that the dust dis- 
tribution properties such as the radial dust distribution or the av- 
erage number of clouds along the equatorial line-of-sight (see 
Honig & Kishimoto 2010, for more details) may change with 
luminosity. This would go in line with the possible trend of 
bluer spectral indices for higher luminosit y objects mentioned 
in Sec t. 15.21 Based on the discussions in iHonig & Kishimoto! 
(2010*), it is possible that higher luminosity AGN either have a 
more compact dust distribution (see also iPoUetta et al.„200&) or 
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have l ess obscuring clouds in the torus (e.g. see lHonig & Becker^ 
Follow-up studies on a much larger sample and luminos- 
ity range should be done to further investigate these trends. We 
will continue this discussion in Sect. 16.31 

A fundamental prediction of all kind of torus models is that 
line-of-sight obscuration, spectral color, and silicate feature are 
correlated. As a first-order approximation, we would expect that 
the SEDs of type 2 AGN are slightly redder than type 1 AGN. If 
we have a face-on view onto the torus, the near-to-mid-lR SED is 
dominated by hot-dust emission from the inner part of the torus. 
In type 2 AGN, the inner hot dust is obscured by cooler dust at 
larger distances from the AGN. This difference in observed dust 
temperature should lead to differences in the spectral color (but 
see Sect. I5.2l i. On the other hand, since hot dust is expected to 
show the silicate feature in emission and cold (optically-thick) 
dust to show a silicate absorption feature, spectral index a and 
silicate feature strength Fs-JFc are supposed to be correlated. In 
the middle panel of Fig. [TT]we plot log Fsi/F^ against the mid- 
IR flux ratio logFv(12.5;um)/Fy(8.5^(m) and spectral index a. 
The type 1 AGN do not show a dependence of the silicate feature 
strength on the spectral index while a range of a = -0.8 . . . -2.2 
is covered. Apparently, the different spectral indices are not an 
effect of obscuration (assuming that the silicate feature depth is 
tracing the line-of-sight obscuration to some degree since the tori 
are seen face-on). Instead it may be a signpost of different dust 
distributions. We showed in Honig & Kishimoto (2010) that the 
spectral slope in the mid-IR is quite sensitive to the power law 
index a of the radial dust distribution //, r" where r is the dis- 
tance from the AGN. The steeper the distribution (i.e. dust more 
concentrated to the inner part of the torus), the bluer the mid-IR 
SED. On the other hand, if the dust distribution is shallower, the 
mid-IR is redder Based on these models, we would conclude 
that the type 1 AGN in our sample show some variety in their ra- 
dial dust distribution, becoming steeper from right to left in this 
plot. We will pick-up this suggestion when modeling our data in 
Sect.|6] 

The type 2 AGN show some weak trend of redder col- 
ors with deeper silicate absorption features. The only excep- 
tion is NGC 7582 were the silicate feature is much deeper for 
the given a than in the rest of the objects. Again, this might 
be connected to additional absorption outside the torus (see 
above). The nominal relation between silicate feature strength 
and spectral steepness for type 2s is logFsi/^c (-0.85 + 
0.20) X logFv(12.5//m)/Fv(8.5//m) (excluding NGC 7582), or 
expressed as spectral index, logFsi/Fc oc (0.14 + 0.03) ■ a 
(Spearman rank p - -0.88, null-hypothesis probability 1.6 x 
10"^). This result is, at least qualitatively, in agreement with 
what we would expect from a clumpy dust torus around the 
AGN. Finally, we can test how well the observed obscuration 
properties in X-ray and the mid-IR are related. In the right panel 
of Fig.[TT] we plot the observed silicate feature strength Fsi/^c 
in the mid-IR against the Hydrogen column densities A^h ob- 
served in the X-rays. Despite some scatter, we find a correlation 
Fsi/Fc oc (0.12 ± 0.03) X logA^H (Spearman rank -0.71, nuU- 
hypothesis probability 9.9 x 10 "*) which is shown as a dashed 
line in Fig. [TT] This is consistent with the analysis of 'S hi et al] 
(l2006l) who used Spitzer data at lower spatial resolution. Fsi/F^ 
can be converted into an observed optical depth TsiUcate in the sil- 
icate feature, and we define Fsi/F^ = exp (-TsUicate) (see right 
y-axis in Fig. [TT] and data in Table |6l). In this way, the correlation 
between observed silicate feature depth and Hydrogen column 
density becomes 

"^silicate — 

(0.00 + 0.73) + (0. 14 + 0.03) x log A^h;22 (2) 



where Nh-ii is A^h in units of 10^^ cm"^. The observed scatter in 
this correlation may be a signpost of a clumpy torus where type 1 
and type 2 AGN are less differentiated. We remind again that the 
observed TsUicate is not the same as the actual optical depth along 
the line-of-sight but reflects the complicated radiative transfer in- 
side the torus. In fact, based on clumpy torus models we would 
expect that the actual optical depth along the line-of-sight is 
much stronger correlated with A^h unless a sizable fraction of the 
Hydrogen column is located inside the sublimation radius (e.g. 
within the broad-line region). Note that NGC 7582 has been ex- 
cluded from this analysis due to the concerns mentioned above. 

Based on the discussion in Sect. 14.31 it is possible that the 
silicate emission features have a different peak wavelength than 
the center of the silicate absorption features. This might poten- 
tially change the outcome of our analysis on the silicate feature 
strength relations. In order to investigate this possibility, we ex- 
tracted F^ JFc at 10.3yum instead of 9.8/vm for the type 1 sub- 
sample. The resulting feature strengths are plotted as gray circles 
in the right panel of Fig.[TT] Changes to the original analysis are 
only moderate and the resulting correlation is shown as a gray- 
dotted line. The results remain consistent within error bars (here, 
Fsi/Fc oc (0.14 + 0.03) x logNn, Speai-man rank p = -0.74, 
null-hypothesis probability 4.5 x 10""*). Thus, our general anal- 
ysis does not suffer significantly from any possible shift of the 
central wavelength in silicate emission features. 

6. Modeling the data with 3D clumpy torus models 

6. 1 . Constraining torus parameters from VISIR N-band 
spectro-piiotometry of nearby AGN 

In lHonig & Kishimotol (l2010h . we prepare interpretation of mid- 
IR observations of AGN using our 3D clumpy torus model 
(Honig et al. 2006). General dependencies of SEDs and inter- 
ferometric visibilities on different model parameters have been 
discussed and it has been shown that the mid-IR SEDs are most 
sensitive to (1) the radial distribution of dust clouds which we 
parametrize as a power-law ?7,.(r ) oc r " (for details please see 
Sect. 2.5 in Honig & Kishimoto 20101) and (2) the obscuration 
inside the torus. The torus-internal obscuration is parametrized 
by the average number of clouds, A^o, along an equatorial line- 
of-sight. We remind the reader that A^o is not equal to the number 
of clouds along the actual line-of-sight unless the torus is seen 
edge-on. We argued that SEDs may serve as a tool to constrain a 
and maybe A^o if the strength of the silicate feature and the spec- 
tral slope in the mid-IR are simultaneously taken into account. 

The mid-IR data presented in this paper has a quite small 
wavelength coverage from approximately 8-13 /vm. This lim- 
its the possibili ties of modeling IR SEDs, e.g. as recently done 
for NGC 1068 (iHonig et alJl2008bl) . As mentioned in the previ- 
ous section, the two main observational parameters that we can 
derive from our VISIR data are the spectral index a of the mid- 
IR continuum and the depth of the silicate feature Fsi/^c (or 
translated into an observed optical depth TsUicate)- In the middle 
panel of Fig.[TT]we plotted both parameters against each other 
and discussed the observed properties of the type 1 and type 2 
sub-samples. Now, we want to compare these observed mid-IR 
properties to our torus models and see if we can constrain pa- 
rameters. For that we simulated a model grid varying a and A^o 
for inclinations from ; = 0° (pole-on) to 90°. As discussed in 
Honig & Kish imotol (1201 0^. other model parameters do not have 
a significant influence on the SED, except for the half-opening 
angle Oq of the torus, which we discuss below. For the moment, 
we use the common assumption of Oq - 45°. 
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Fig. 12. Comparison of observed mid-IR properties of our type 1 and type 2 AGN to model SEDs simulated with our 3d clumpy 
torus model. Left: F^i/F^ and TsiUcate of the type 1 sub-sample (blue triangles) plotted against the spectral spectral index a of the 
mid-lR continuum emission. Overplotted are model predictions for parameters ranging from a - 0.0 (orange) to a = -2.0 (light 
blue) and A^o =2.5 (asterisks) to A^o = 10 (x-shapes) for a mean type 1 inclination of / - 30°. The black arrows note the directions in 
which the model points change when varying a. No, and /. Right: F^-JFc and TsiUcate of the type 2 sub-sample (red squares) plotted 
against the spectral spectral index a of the mid-lR continuum emission. Oveiplotted are model predictions for parameters ranging 
from a - 0.0 (orange) to a = -2.0 (light blue) and A^o - 2.5 (asterisks) to A^o - 10 (x-shapes) for a mean type 2 inclination of 
i - 75° (colored symbols) and / - 60° (gray symbols, illustrating the inclination effect). The gray arrows are a rough illustration as 
a guide for the reader in which direction the model points m ostly change when i ncreasing a, Nq, and /. For reference, NGC 1068 is 
shown as an orange star, based on the mid-lR spectrum from lMason et al.l (l2006l) . 



In Fig. [12] we show model predictions for the silicate feature 
strength Fs-JF^ (or observed optical depth TsiUcate in the silicate 
feature) and the spectral index a of the mid-IR continuum. In 
the left panel, we show simulations for an inclination angle of 
; - 30° which reflects a typical type 1 AGN line-of-sight for a 
half opening angle of 6o = 45°. Model results for different ra- 
dial power law indices are color coded from a - -2.0 (blue) to 
a = 0.0 (orange). For each a, several values for A^o have been 
calculated and are marked by different symbols in the plot from 
A^o = 2.5 (asterisks) to A^o = 10 (x-shapes). Dashed-colored lines 
connect symbols with the same a but varying A^o, while dotted- 
black lines represent the same A^o for different a values. The right 
panel of Fig.[T2]shows the same model grid but for a typical type 
2 inclination of / = 75°. We also overlaid the grid for / = 60° to 
illustrate how inclination effects change the results. In general, 
changing a has a strong effect on t he spectral index a (as dis- 
cussed in iHonig & Kishimotoll201(]h and some effect on Fsi/F^. 
On the other hand, varying A^o and / has almost no effect on a 
but changes F^JFc- It is noteworthy that the spectral index a 
is very insensitive to the torus inclination, even when chang- 
ing from edge-on to face-on line-of-sights. As discussed, this 
is characteristic for clumpy models. The gray arrows in Fig.fT2l 
are intended to help the reader, summarizing the rough trends of 
how the different parameters change the place of an object in the 
model grid. 

The model grids in Fig. [12] illustrate that some fundamental 
parameters of the dust torus can be constrained even using the 



limited wavelength range covered by our observations. We over- 
plot the type 1 and type 2 sub-samples in the left and right panel 
of Fig. [12] respectively, in order to see what region in parameter 
space is occupied. As already discussed, the most robust con- 
straint can be put on a since it is the almost exclusive parameter 
to determine an object's place in horizontal direction. The type 
1 AGN in the left panel populate the range between a - 0.0 to 
a = -1.5, with some clumping roughly at -1.0 + 0.5. ESO 323- 
G77 is the only exception being located between a = -1.5 and 
a = -2.0 but outside the displayed model grid. It is discussed in 
more detail below. Interestingly, most type Is also cover a rather 
naiTow range in vertical direction, somewhere between No - 5 
and 7.5, and they follow approximately a line of equal A^o for 
varying a. NGC 7213 is slightly above the other objects, sug- 
gesting a lower value for A^o- Since it is not a classical Seyfert 
galaxy but a LINER, the offset from the Seyfert Is might indi- 
cate some intrinsic difference in the dust structure. On the other 
hand, it may also result from an inclination effect: If NGC 7213 
is seen more pole on than the other objects, then we would expect 
that for the same a and A^o the object is placed slightly above the 
others (see also inclination effect illustrated in the right panel of 
Fig.fTSTi. However, at low inclinations, any viewing-angle effect 
is small and not able to explain all of the offset of NGC 7213. 

The right panel of Fig. [12] shows the type 2 AGN of our 
sample. Most objects are again located between a - 0.0 and 
a = - 1 .5 and they follow the line of constant A^o - 5 for in- 
clination / - 75° or A^o ~ 6 for / - 60°. There are, however, a 
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number of outliers from this plot, most notably NGC 7582. For 
this galaxy we already remarked that additional extinction from 
a circumnuclear star-burst ring might interfere with the emission 
from the nucleus. Consequently, the observed mid-IR proper- 
ties are not generic properties of the torus emission. To study 
the possibility if such additional (screen) absorber may be re- 
sponsible for outlier in the plot, we used the extinction curves 
of our dust composition (I SM dust with Ossenkopf silicates; see 
iHoni^ & Kishimotol2010l for details) and derived the change of 
mid-IR spectral properties for Ay = 10. The resulting shift is 
shown as an orange arrow in both panels of Fig. [121 Actually 
NGC 7582 and MCG-5-23-16 - the objects in our sample 
which are located below the model grid - have a comparably 
high inclination of the host-galaxy which may contribute some 
extinction screen to the nuclear emission. On the other hand, 
MCG-5-23-16 is very close to the model grid and the offset can 
be easily explained by either a torus inclination / > 75° or a 
torus with A^o - 10 or slightly higher. Please note that "redden- 
ing" (extinction) is actually making the A^-band mid-IR colors 
slightly bluer, which is an effect of the extinction curve at these 
wavelengths. Another interesting object is NGC 5643. This low- 
luminosity AGN has the reddest colors in our sample. It follows 
quite well the line of constant A^o occupied by the other AGN, 
but is placed at a > 0.0. This would imply that the torus is much 
more extended than in the other objects, compared to the intrin- 
sic scale, something that might be tested by IR interferometry. 

Finally, we will discuss ESO 323-G77 which is an outlier in 
the type 1 plot (left panel of Fig.fTSli. It is located slightly below 
the model grid between a = -1.5 and -2.0. Based on our mod- 
eling, the offset cannot be explained by either an inclination ef- 
fect or a higher N[). Inclination is excluded since ESO 323-G77 
is even below the type 2 AGN model grid in the left panel of 
Fig. [12] Moreover, for such steep dust distributions, most clouds 
are located in the innermost part of the torus close to the sub- 
limation radius. Adding more clouds to the torus will dramati- 
cally increase the torus-internal obscuration and moves the mod- 
els rather to the right than downward in the plot. This effect can 
already be seen in the model grid for a - -2.0 and A^o = 10 
and becomes even stronger for higher A^o- On the other hand, 
ESO 323-G77 is showing some absorption pattern in the optical 
and the X-rays (see Table [Tj) which might result in a down-left 
drift in the plot (see orange arrow). It is located in the same area 
of the Lmir - Lz-correlation as MCG-3-34-64 and IC 5063 for 
which we found extended emission probably coming from the 
NLR. Since ESO 323-G77 is a type 1 AGN with a more or less 
face-on line-of-sight, dust in the NLR can cause screen extinc- 
tion which may explain the offset position in Fig. [121 In this 
case it is reasonable to conclude from Fig. [12] that the intrin- 
sic dust distribution power-law index of ESO 323-G77 is most 
likely somewhere between a - -1.5 and -2.0 as expected from 
the relatively blue mid-IR color 

In summary, based on our modeling we find that the Seyfert 
galaxies in our sample have a torus with a rather shallow ra- 
dial dust distribution with a typical power law index of approxi- 
mately a - -1 + 0.5. This parameter is quite solidly constrained 
by the mid-IR spectral indices. For the presented model grid, we 
also find a typical range for the average number of clouds along 
an equatorial line-of-sight of A^o ~ 5-8 for the whole sam- 
ple. It is, however, difficult to pin down this parameter for each 
individual object since inclination effects can slightly alter the 
results. As of yet, we did not discuss how the vertical structure 
of the torus (represented by the half-opening angle, Oo, of the 
torus) influences the models. In fact, does not have an effect 
on the horizontal position of the object (exclusively determined 
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Fig. 13. Comparison of the VLTI/MIDI 8 - 13/zm interferom- 
etry of NGC 3783 with our 3D clumpy torus models. The in- 
terferometr y data at projected b aseline length 65 m have been 
repo rted by Beckert et al.l (l2008h . the 43 m data was presented in 
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(I2009ah . We plot the observed 8-13 jim visibil- 



ity (color-coded from red to blue) aga inst the spatial wavelength 
in units of the inner radius (see Ki shimoto el al.i 2009a, for 
details). Overplotted are the model visibilities for 8.5 yum (red 
line) and 12.5 fj.m (blue line) based on model parameters derived 
from Fig.[T2](fl = -0.75, A^o = 7). The 8.5 and 12.5 /im data 
have the same colors as the corresponding model. 



by a) but can change its vertical position, i.e. the strength of the 
silicate feature. We find that increasing Oq from 45° to 60° has a 
similar effect as decreasing the inclination angle by about 15°. If 
we consider 9o - 45° as a reasonable sample average, the range 
of A^o ~ 5 - 8 found as a typical value for the whole sample is 
still valid, although it might be complicated to constrain A^o, 0q, 
and i for each individual object. 

6.2. Comparison to IR interferometry 

While our modeling approach based on the spectral slope and sil- 
icate feature strength is very limited in wavelength range, we can 
test its results and predictions. For that, we compare the derived 
radial dust distribution to recent results from IR interferometry. 
In Kishimoto el al. ( 2009a), we argued that IR interferometry of 
type 1 AGN can be used to determine the surface brightness pro- 
file of the dust torus. For that, the A^-band visibilities observed 
by VLTI/MIDI have been plotted against an object-intrinsic spa- 
tial scale (the near-IR reverberation radius, which is presumably 
indicative for the inner boundary of the dust distribution). It has 
been found that the brightness profile in NGC 3783 depends on 
the radius oc r"^. We can now try to reproduce this finding by our 
models using model parameters as constrained by Fig. [12] For 
NGC 3783 we find that a ~ -0.75 and A^o ~ 7 for / = 30. Based 
on these parameters, we simulated images and calculated corre- 
sponding mid-IR visibilities. The results at 8.5 fim and 12.5 fim 
are shown in Fig. [13] The model visibilities are in general agree- 
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ment with the available data, despite some overestimation of the 
visibility at longer wavelengths for the long projected baseline 
(baseline length 65 m, PA 120°). The fact that the observed vis- 
ibilities are lower in PA 120° may be due to some elongation 
of the source. Based on our clumpy torus models we would 
not expect that the torus itself produces these kind of deviation. 
Interestingly the larger size in PA 120° is pointing into a simi- 
lar di rection as the polarization angle (0poi ~ 135° Smith et al. 
I2OO2I) which implies that the mid-IR nucleus is elongated to- 
wards the NLR. This may place NGC 3783 into the same group 
of objects as IC 5063 and MCG-3-34-64 which show some con- 
tribution from the NLR to the mid-IR emission (although for 
NGC 3783 on much smaller scale and fraction), in particular to 
longer wavelengths. We leave a more detailed analysis of the 
interferometric measurements to a much larger dataset that was 
obtained recently (Honig et al. in preparation). Despite these de- 
tails, NGC 3783 as a case study illustrates that the brightness 
profiles of the dust distribution are strongly connected with the 
mid-IR slope a, and both are very sensitive to the radial distri- 
bution of the dust. 

In the left panel of Fig. [121 we also show NGC 1068 as an 
orange asterisk. It is located slightly below the / = 75° model 
grid in the range of a = -1.0 ... - 1.5 and A^o ~ 10. Recently, 
we presented detailed modeling of the full IR SED of NGC 1068 
at high-spatial resolution simultaneously with near- and mid-IR 
interferometry (Honig etal. 2006, 2007, 2008b). Since IR in- 
terferometry resolves the spatial brightness profile, it is consid- 
ered as very constraining for model parameters. For NGC 1068 
we found typical parameters of a ~ -1.5, / = 70 - 90°, and 
A^o ~ 10 which is in excellent agreement with the analysis pre- 
sented here based on the mid-IR spectral index and the silicate 
feature strength. 

In summary we find that the results from our analysis of the 
mid-IR spectral index a and the silicate feature strength is in 
good agreement with results from IR interferometry. Model pa- 
rameters derived by our spectro-photometry approach are con- 
sistent with a much more detailed analysis of the IR emission of 
NGC 1068 and can reproduce interferometric measurements of 
NGC 3783. In turn, based on the analysis of the mid-IR spectra, 
we can make predictions about the spatial brightness profiles of 
AGN for interferometry. For instance, we expect that ESO 323- 
G77 and MARK 509 have quite compact profiles (i.e. high vis- 
ibility) while objects Uke NGC 3227, NGC 7213, or NGC 7469 
have rather extended profiles (i.e. low visibilities) when scaled 
for their intrinsic radii. 

6.3. A luminosity dependence of tiie torus properties? 

One interesting aspect of the torus is whether or not there are 
any differences in its structure at low and high AGN luminos- 
ity. Obscuration statistics of type I and type 2 AGN suggest 
that the r elative fraction of type I AGN increa s es with lu minos- 
ity (e.g. Simpson 20051: iMaiohno et al.ll2007L iHasineed [20081) . 
This phenomenon is commonly explained by a decrease of torus 
covering factor with luminosity - and referred to as the "reced- 
ing torus". On the other hand, if the covering factor decreases, 
the "reprocessing ratio" Li2f,m/Lx ratio is also expected to de- 
crease since less torus surface is heate d by the AGN . How ever, 
based on the Lmir - Lx -correlation bv lGandhi et al.l (1200 9). we 
find that L^^mlLx °^ jj^'^^^oi which means that the reprocess- 
ing ratio does not significantly change with luminosity (maybe 
marginally increases; but we note that when using the bolometric 
luminosity instead of Lx, the ratio might slightly decrease with 
luminosity, depending on the bolometric correction used; see 



also lGandhietan(l2009t) . Sect. 4.5). In lHonig & Beckerj (l2007h . 
we suggested a possibility to explain these observations: increas- 
ing AGN radiation pressure may cause large dust clouds to be 
driven out of the torus at higher luminosities. Consequently, the 
absorbing column towards the AGN decreases and it becomes 
less likely that the line-of-sight towards the nucleus is obscured. 
On the other hand, the covering factor and reprocessing ratio will 
not change. 

The question remains, however, if there are direct hints that 
the torus dust distribution changes with luminosity. The most 
direct way of probing this possibility is, of course, IR interfer- 
ometry. On the other hand, due to the sensitivity of current in- 
terferometers, one is limited to few objects and the luminosity 
coverage is limited. In the previous sections we showed that the 
radial structure of the torus and the torus-internal obscuration 
properties can be constrained by simultaneously comparing the 
mid-IR spectral index a and the silicate feature strength. In par- 
ticular, constraining the radial power-law index a using a seems 
to be quite solid (see also Sect. 16.2b . 

The AGN in our sample have typical luminosities of nearby 
Seyfert galaxies, i.e. in the range of 10"^^erg/s to 10"*^erg/s, 
and conclusions drawn about the typical torus properties of 
a - -I + 0.5 and A^p - 5-8 are on ly valid for these lu- 
minosities. Recently, iPolletta et al.l (l2008h presented rest-frame 
near-to-mid-IR SEDs of 23 type 2 QSOs with luminosities be- 
tween lO^'^erg/s to several lO'^^erg/s. Interestingly, most of 
them showed quite blue IR spectral indices. Using the clumpy 
torus models, it was concluded that the typical radial power-law 
indices in the type 2 QSOs were tracing a quite steep dust dis- 
tribution, with a ~ -2 ... - 3. If this is indeed a typical value 
for luminous AGN, it is much steeper than what we find in our 
intermediate-luminosity sample. Consequently, there could be a 
change of the torus structure with luminosity. 

At the low-luminosity end, we do not have a good cover- 
age of mid-IR data where spatial resolution is sufficient to iso- 
late the torus emission. However, we can compare the lower- 
luminosity objects in our sample to those at the higher lumi- 
nosity end. In fact, in Sect. 15.21 we note that the reddest ob- 
jects in both type 1 and type 2 sub-samples are those with the 
lowest luminosities (NGC 3227 and NGC 7213 for type Is, 
NGC 5643 and ESO 428-G14 for the type 2s). This would sug- 
gest that the lower-luminosity AGN have much more shallow 
dust distributions, i.e. the torus is relatively more extended as 
compared to higher-luminosity objects. If true, we might see an 
imprint of this behavior in the mid-IR size-luminosity relation, 
e.g. as it can be constructed from interferometry. The prediction 
would be that for a sufficiently large sample the relation does not 
scale as r oc L}I^ but with a slightly smaller power-law index. 
Similarly, near-IR interferometric sizes would be systematically 
larger than IR reverberation-mapping radii at lower luminosi- 
ties, and would better matc h each other at higher luminosities 
(cf. Kish imoto el al.ll2"009bl for Keck interferometry results of 4 
objects). 

7. Summary 

In this paper, we presented high-spatial resolution mid-IR 
spectro-photometry of a sample of nearby AGN. Based on the 
selection criteria, the sample cannot be considered complete but 
the objects are typical for local intermediate-luminous AGN, 
i.e. Seyfert 1 and Seyfert 2 galaxies ranging from Compton-thin 
to Compton-thick X-ray obscuration. The ground-based mid-IR 
data was taken with the VLT mid-IR imager and spectrograph 
VISIR, and it represents, to our knowledge, the largest sample 
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of mid-IR A^-band spectra of AGN with angular resolution <0'.'4. 
The main characteristics of these data are: 

• Our ground-based observations isolate the nuclear dust emis- 
sion of the torus from host-galactic sources. PAH emission 
features, commonly associated with star-formation in the 
host galaxy, are mostly absent in our data (i.e. at spatial 
scales <100pc), although seen in Spitzer data at lower spatial 
resolution. We showed that subtracting the nuclear VISIR 
emission from the Spitzer data reveals spectra reminiscent of 
a star-formation template. 

• The mid-IR spectra show only moderate silicate features. 
While silicate absorption is easily identified in type 2 AGN, 
the corresponding emission feature in type 1 AGN is very 
weak, if present at all. Thus, the weakness of the silicate fea- 
ture is an intrinsic property of the emitting region which is 
<25-160pc, depending on the object. This points to an ori- 
gin in the circumnuclear dust distribution - the "dust torus". 

• The strongest silicate emission features are seen in the type 2 
AGN NGC 2110 (see also Mason et al. 2009) and the broad- 
line LINER NGC 7213. We analyzed the silicate emission 
feature of the latter object and found a shift of central wave- 
length of the silicate feature to ~10.5 yum. This shift is prob- 
ably not exclusively caused by radiative transfer effects due 
to absorption as recently suggested, but might indicate dif- 
ferent dust composition than in the standard ISM, or a radial 
change of grain size and/or chemistry in the dust distribution. 

The observed mid-IR luminosities have been compared to 
other AGN luminosity tracers, namely the 2-lOkeV luminosity 
Lx and the luminosity of the [O iii](/l5007 A) narrow emission 
line. We showed that correlations of both tracers with Ln^m are 
quite strong and tighter than between Lx and L[o m] . It may be 
worthwhile considering Ln^m as an isotropic AGN luminosity 
tracer. 

Three object in our sample, MCG-3-34-64, ESQ 323-G77, 
and IC 50 63, are outliers in the otherwise tig ht Lmir - Lx- 
correlation (iHorst et alJl20()8tlGandhi etani2009h . We evaluated 
mid-IR images in 4 different filters per object and found that they 
show significant extension and elongation. In IC 5063, the elon- 
gation direction is consistent with HST [Oiii] emission maps. 
We suggest that a sizable fraction of the mid-IR emission in 
these objects is originating in the narrow-line region, possibly 
caused by dust located there, similar to the mid-IR-bright NLR 
in NGC 1068. This additional emission may be the reason for 
the mid-IR "overluminosity" in the Lmir - Lx-correlation. 

Since the wavelength range covered by our data is limited 
to 8 - 13 fim, the main properties that we can use for analyzing 
the dust emission are (1) the mid-IR continuum spectral index a, 
and (2) the strength of the silicate feature. We found that there 
are at best mild trends of increasing a with increasing Li2/jm and 
increasing a with increasing Hydrogen column density A^h- The 
weak trend of the a vs. A^h shows that strong or weak obscura- 
tion does not have a significant impact on the overall color of the 
mid-IR SED, consistent with the as yet not detected difference 
bet ween the Lmir -Lx-co rrelation of type 1 and type 2 AGN (see 
also Gandhi et a n 12009^. As suggested in Honig & Kis himota 
(1201 OK the distribution of the dust around the AGN has much 
more impact on the mid-IR spectral slope than obscuration or 
orientation effects. We also compared the strength of the silicate 
feature f si/^c (or observed TsiUcate) to the mid-IR luminosity and 
the Hydrogen column density. We find a similar correlation of 
silicate feature strength and A^h as Shi et al. (2006), correspond- 
ing to Tsiiicate 0.14 X logA^H- Interestingly, there might be a 
difference between silicate feature strength at lower and higher 



luminosity in our sample: type 1 and type 2 objects at lower 
Li2/jm seem to show a stronger difference in silicate features than 
objects at higher Ln^m where type 1 and type 2 features appear 
rather similar. A larger sample of objects with high angular res- 
olution is required to see if the trend can be confirmed. 

Finally, we plot silicate feature strength against mid-IR spec- 
tral index a. Such an analysis is motivated by our 3D clumpy 
torus model where we found that simultaneously accounting for 
silicate feature and spectral slope may be very constraining for 
torus parameters. We overplot the observed properties of type 1 
and type 2 AGN with predictions from our model, varying (1) 
the radial dust distribution power law index a (radial distribution 
rjr oc r°), (2) the average number of clouds A^o along the line-of- 
sight in the equatorial plane, and (3) the torus inclination Our 
main conclusions are: 

• a versus silicate-feature-strength plots are quite constraining 
for torus parameters since a is almost exclusively sensitive 
to a and the silicate strength is mostly sensitive to A^o and 
/. We discussed how the position of an object is depending 
on these and possible other model parameters and concluded 
that a can be constrained for each individual object while A^o 
may be better derived as a sample mean. 

• For our Seyfert sample, we find typical values of a = -1.0 + 
0.5 and A^o = 5-8. The sample average values are similar for 
type 1 and type 2 AGN, which is a support for the unification 
scheme. We note, however, that this conclusion is only valid 
for typical Seyfert AGN as presented in our study. 

• Some objects may be suffering from extinction in the host 
galaxy. In such cases, however, additional screen absorption 
changes mostly the object's vertical position in the a-silicate 
strength plot. 

• The approach of constraining a and A^o based on the mid- 
IR spectral index a and the silicate features strength has 
been tested for examples where we also have much broader 
IR SEDs and/or infrared interferometry (NGC 1068 and 
NGC 3783). Interferometry directly traces the brightness dis- 
tribution of the emission region and is most sensitive to the 
radial distribution of the dust (see Kishimoto el al. 2009al). 
We find that the model parameters derived from simultane- 
ous SED and interferometry modeling of NGC 1068 are con- 
sistent with parameters derived from the a-silicate feature 
plot. In addition, a and silicate-feature-strength parameters 
derived for NGC 3783 are able to reproduce the available 
mid-IR interferometric data. Based on our mid-IR single- 
telescope data, it is possible to make predictions for inter- 
ferometry. 

Finally, we discuss a possible change of torus properties with 
luminosity. Based on a comparison with literature data and mod- 
eling, we find evidence that the dust tori might be more compact 
(i.e. a X -2 ... - 3) for high luminosity AGN and shallower 
on the low-luminosity end (a w 0.0 ... - 0.5), but we need bet- 
ter selected samples to confirm this trend. If true, this change 
of the dust distribution in the torus will slightly alter the mid- 
IR size-luminosity relation and causing differences when com- 
paring near-IR interferometric sizes with near-IR reverberation 
mapping radii of AGN. 
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